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Poiseuille Flow of Certain Anisotropic Fluids 


J. L. ERICKSEN 


1. Introduction 

Simple shearing flow of relatively simple, incompressible anisotropic fluids 
is analyzed by ERICKSEN [J]. For a certain class of these, the predicted behavior 
is like that of the so-called Bingham materials. That is, the steady state shear 
stress is a linear but not homogeneous function of the rate of shear. The purpose 
of this paper is to present analyses indicating how these might behave in Poiseuille 
flow, t.e. in steady flow through an infinitely long circular pipe. It turns out 
that there are infinitely many solutions which seem more or less acceptable. 
Acceptance of a plausibility argument rules out all but one on the grounds 
that the remainder represent unstable situations. That which remains yields 
results entirely similar to those deduced by BucKINGHAM [2] and REINER [3] 
for Poiseuille flow of Bingham materials. Though this might be expected, it is 
not entirely obvious a priort. In the one dimensional descriptions of Bingham 
materials used by BUCKINGHAM and REINER, or in the three dimensional formu- 
lation introduced by OLpDRoyD [4], two sets of constitutive equations are used. 
One set, applying at low stresses, describes an ordinary elastic material. The 
other, applying at high stresses, involves a viscosity similar to that which 
characterizes Newtonian fluids. A scalar yield criterion is used as a criterion 
for determining which set applies. The theory used here employs a single set 
of constitutive equations and no yield criterion. One consequence is that the 
reasoning used in analyzing Poiseuille flow of Bingham materials is not directly 
applicable to the materials here considered. 


2. Equations of Motion 
It is a straightforward matter to refer the equations proposed by ERICKSEN 
[1, eqs. 19 and 20] to cylindrical coordinates (7, #, z), then to introduce the 
assumptions that the velocity field is of the form 


dy Cie) ihe 
(2.1) = ap oe ge I 


and that all quantities are independent of #. One obtains the stress-deformation 
relations EBON PEST nn, 
leg =— Pp+X 5 
t.a=—p+X nit 2usf n,n, 
(2.2) t,, =X n,n, + [ug(t — 5) +] 7’, 
tp = XN, Ng + Ugh Np MN, 
tg, = X NgN, + Ugh M, Np, 
X = py t+ bef n, M,. 
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Here t,,, ty», etc., are physical components of the symmetric stress tensor, p is 
an arbitrary pressure and (n,, %,”,) are physical components of a unit vector 
pointing out a preferred direction in the fluid. Greek letters denote material 
constants. The vector satisfies the differential equations 


dl 
(2.3) Gnd — — Afi ny ng M,, 
2 Be = f'n, (A+1—2An?), 
where : 
=F ! re 


is the material derivative operator. The fluid behaves as a Bingham material 
provided 

(2.4) [Atel s jae 0; 

which we assume. 


Flows of the form (2.1) are accelerationless. In the absence of body forces, 
the stresses must thus satisfy the equations of static equilibrium, which reduce to 


et ae a 
pe ae 4 ; (ty — toa) =90, 
Ct, 9 Ct 2 
(2.5) je vi vi bia Os 
ot Nie 1 
oy a 1 Vv ba =O 


From (2.1), df'/dt=0, so /’ may be treated as a constant in integrating (2.3). 
Equations of precisely the same form are analyzed by EricKsEN [J, 5], from 
which we may conclude that, when (2.4) holds, essentially all solutions of (2.3) 
approach the solution 


(2.6) m=SNQ, M=0, n,=cosp 


as too. Here is a constant angle determined by the conditions that 


(277) tan?» = (A — 1)/(A+1), 

and that it lie in the range 

(2.8) 0<p<a/4 WROD: gfe Oss amends 
(2.9) 32/4<y<a when f’<0, A>1, 
(2.10) mf A< gia] 2h  whenay —<20,. 9 Mies ee 
(2.11) R2< P< 3a/4 when) 705 ft 


The theory treats m and —vn as indistinguishable, which means that there is 
no loss in generality in taking 0< @<z, as is done here. 
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; 3. Steady State Flow 


We use the stable, steady state solution of (2.3) given by (2.6)—(2.11). 
i proceeding under the assumption that /’=-0. Using (2.2), we then find that 


(3.1) t,=—P+msinPo+Af, 
(3.2) tgg=—?, 

, (3.3) t,,=—ptmco?o+Bf, 
(3.4) t,,=C+Df, 

| (3.5) to =te,=0, 


) where A, B, C and D are constants given by 


(3.6) 2A = (My sin® p + 2y,) sin2, 
(3.7) 2B = (ul, cos? y + 2g) Sin 29, 
(3.8) 2C— 7, sin. 20, 

(3.9) 4D = 4 (lg + ) + fy Sin? 29. 


Since ¢,,>C as f’->0, C is interpretable as an apparent yield shear stress. 
Similarly, D is an apparent fluidity. It should be noted that the constants 
A, B and C reverse sign when /’ does, so places where /’=0 must be treated 
with care. For the moment, we exclude these. It is then a simple matter to 
integrate (2.5) to obtain 


(3.10) p=msin'g + Af’ + {+ (14 sin?» + Af’) dr — Ez, 
(3.11) ty,=C+Df' =—Er2+Fr, 


where £ and F are constants of integration. 
Physically, we expect the stresses to do non-negative work in deforming 
the fluid, 


(3.12) t,f =(C+Df)f2o. 
From (3.8), (3.9) and (2.7)—(2.11), we find that this requires that 


(3.13) eee, fg (A2— 1) + 402 (ug +m) 20 when A>1, 


(3.14) MyS0, fp (A2—1) + 4A?(ug +m) 20 when A<—1. 


We assume that the strict inequalities hold in (3.12) and (3.13). 

It is easily seen that our assumptions rule out the possibility of having a 
solution devoid of singularities at r=0. Previously mentioned analyses of 
Bingham materials suggest that, near the center of a pipe, plug flow occurs. 
That is, a central core moves as a rigid body. We now explore this possibility. 


4. Plug Flow 


We seek a steady solution of the basic equations with f/=/)—const. which 
is devoid of singularities at y=0. From (2.3), the material derivative of nm must 


4* 
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vanish, which implies that n=n/(r, 8, z—fot). Dependence on # was excluded 
earlier and, except when f,=0, dependence on z— fq ¢ is ruled out by the assump- 
tion of steadiness. We thus assume that 


(4.1) nN, =N,(7), TG My (7) , Ue n,(7). 
From (2.2) and (2.5), we then have 
rt, = 4 0, Mg = G = const, 


For ¢,, to remain finite as y->0, G=O, so that 


(4.2) N, Ny = 0. 
From (2.2), 
(4.3) ip — ON 2 Oe 


As will be discussed later, (4.3) makes the possibility ,=0 uninteresting for 
pipe flow. Since m is a unit vector, (4.1) and (4.2) leave us with the possibly 
interesting alternative 


(4.4) nh, = SI, <g= 0, B= cosy, 


where yw is any function of 7. From (2.2), we then have 


(4.5) t,=—p+msin?y, 
(4.6) isg7= — P, 

(4.7) t,,=—pt+mcos?y, 
(4.8) t,,=4m,sin2y, 

(4.9) t,o = ty, = 0: 


Integrating (2.5) subject to the requirement that ¢,, remain finite as y->0, we 
obtain 


(4.10) pP=msin?y + iy [ sin*pdr — Hz, 

(4.11) 2t,,=msin?'y = — Hr, 

where H is a constant of integration. From (4.9), we should have yO as 
v—>0 to keep the normal stresses finite as y—>0. Because of the indistinguishability 


of n and —n, there is no loss in generality in taking OSy<a. Solving (4.11) 
for sin py subject to these conditions, we obtain 


(4.12) sin?y =} [1 — V1 — (A7/)?], 
with 

(4.13) OSysa7/4 if Hiu<o, 
(4.14) 3nAS ps it Bla 0; 


This gives solutions of the desired type, valid for 
(4.15) OS 77 =| /A\. 
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5. Combined Flow 

We consider the possibility of using a solution of the type discussed in §3 
for r2r,, a solution of the type discussed in §4 for y<7,, requiring that the 
velocity vectors and stress vectors should match at r=7,. That is, f, t,,, tye 
and #,, should be continuous at y=7,. For the moment, 7, represents any value 
j of y in the range O<7X~,, where 7, satisfies (4.15). From (3.5) and (4.9), t,» is 
| automatically continuous, being zero according to both solutions. 

For plug flow, we have from (4.5), (4.10) and (4.12) 


tL 


my tr) = a eae es (H xjuy)*|dx + Ha + K, r<1, 


0 
where K is a constant of integration. From (3.1) and (3.10), we may write 


(5.2) tn) =— [ L(msintp + Af)dx+Ez+L, r= 1p. 
For (5.1) and (5.2) to agree for r=y7,, z arbitrary, we must have 
(5.3) Es; 


Yo 


(5.4) L=—F yy [ -[—Vi— Hxfuy)*] dx +k. 


0 


Using (5.3), (3.11) and (4.11), we see that ¢,, is continuous at y=7, if and only if 


(5.5) F=0 
which means that 
(5.6) t,,= — Hr/2 


throughout. From (3.11) and (5.6), 
aes oe 2 
(5.7) eae eh LAC?) M, 
where M is an arbitrary constant. For f to be continuous throughout, 
(5.8) f= — ap (H+ 4Cn) + M. 


As long as we satisfy these conditions, we obtain a combination of solutions 
which is physically acceptable for all values of 7. This would leave us with an 
overabundance of solutions to the Poiseuille flow problem. We introduce a 
plausibility argument which rules out most of the possibilities on the grounds 
that they describe unstable situations. 

From the plug flow solution, we can calculate ”, the limiting value of n 
for 7-7, 7<Yg, as a function of 75, 


(5.9) n,=siny (7), %=0, %,=Cosp(re). 
For the flow external to 7,, we have from (5.7) the limiting velocity gradient 


(5.40) f(t) =— sp (An +20), 
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which will in general tend to rearrange the structure represented by n. That 
is, if we regard ” as an initial value of m in (2.3), we see that f’(72) would cause 
n to change in accordance with 


2 Gtr = f(y.) n,(A —1— 2An?), 
ny as = fi (ro) n,(A+1—2An}3), 


nm» remaining zero. We expect this to make the interface unstable in general. 
Exceptions might occur if 


(5.41) Y (to) = 9, 


when n coincides with the stable, steady state value of n in the flow exterior 
to 7=75, oF if 


(5.12) f'(r2) = 0 
which implies that dn/dt=0. If (5.11) holds, (4.11) and (3.8) yield 
(5.13) 2 = — (t4/H) sin 2y(r.) = — (m/H) sin 29 = — 2C/H, 


whence it follows from (5.10) that (5.12) holds. Hence (5.12) must hold and, at 
the interface f SSL ee 
Thus the shear stress takes on its apparent yield value at the interface. This 
is in accord with analyses of Poiseuille flow in Bingham materials, wherein it 
is assumed or deduced from other assumptions that, at the boundary of the 
rigidly moving plug, the velocity gradient vanishes, the shear stress taking on 
its yield value. 
If (5.12) holds, (5.13) must hold, from which we may infer that either (5.11) 
holds or 
(5.14) p (72) = 2/2 — @. 
Noting (4.13) and (4.14), then referring to (2.8)—(2.11), we see that (5.11) but 
not (5.13) can hold if A>1, whereas (5.13) but not (5.11) can hold if A<—1. 
The rather crude argument used here hardly settles the stability question 
for these flows. Among other things, it says nothing concerning the stability 
of the noted exceptional cases. It seems probable that, under some conditions, 
this exception describes a stable situation, but we leave this question to the future. 


6. Poiseuille Flow 


We consider using the solutions developed to describe steady flow through 
an infinitely long, stationary circular pipe of radius R, to which the fluid adheres, 
regarding as given the pressure head 


H = — 0p/0z = dt,,/dz. 


Assume the consequences of the plausibility argument used in §5 and suppose 
H is sufficiently large to make 


(6.1) ¥,=—2C/H<R. 
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Using the adherence condition to evaluate the constant M in (5.7) and (5.8), 
we have 


f=—- Zr —R) [A +R)+4C], 21, 


f=f=~yle— RH +R)+4C], Sn, 
where C and D are given by (3.8) and (3.9), viz 
2C=m,sin2g, 4D=p.sin*29+4(us+y). 
To determine the appropriate value of y, we note that (5.13) requires that 
—mHsin2y>0, 


Then use (3.13) or (3.14) and (2.7) —(2.11). For a given material, we thus uniquely 
determine the velocity distribution for the combined flow as long as (6.1) holds. 
With the understanding that C and D are interpretable as yield stress and 
fluidity, it is indistinguishable from the velocity distribution for a Bingham 
material. We thus obtain the same expression for mass flux in terms of pressure 
head, given e.g. by REINER [3]. 

If the pressure head is sufficiently low, (6.11) must fail, the radius 7, then 
lying outside the pipe. A reasonable interpretation is that the material remains 
at rest. In other words, the plug flow solution with f,=0 applies. From (4.5) to 
(4.14) and from (5.1) we then obtain the stress 


Pagan ae 
tog = Y —Fuy(1 \1- (H y/u)"], 


(6.2) tey =Y + My | Ase (H r/y)?, 


where ‘ 
v= sin f [4 Vi — (Hxfy)*] dx + Hz +K, 
0 
when flow occurs, these stresses obtain for O<rS7,. For 7,5 7< R, we then 
similarly obtain 


t= Dis 
typ =Z — py sin? + (Hr +20), 
A—B 
(6.3) tag =Z Tt ty COS 29 + (74-2), 
t.,=—+Hr, 
t,9 ='g,=0 


where A, B, C, D are constants given by (3.6)—(3.9), 


(6.4) Z=— [ L(msin'g + Af)dx+Ez+L, 


' 
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and L is given by (5.4). Carrying out the integration in (6.4), we obtain 
Z = (ACID — yy sin® g) In(o/r,) +S 7 — me) +EZ+L. 


Except for K, which remains arbitrary, all constants are uniquely determined 
as functions of H. Changing K changes the stress by a uniform hydrostatic 
pressure, which has no effect on incompressible materials. 

The stresses given by (6.2) or (6.3) differ from those predicted by OLDROYD’s [4] 
equations for Bingham materials, which imply that t,,=ts,—t,,. However, in 
terms of quantities generally discussed in connection with Poiseuille flow of 
Bingham materials, the two theories yield indistinguishable predictions. 

For Poiseuille flow, the possibility indicated in (4.3), 7.2. n,=0, %»== 0, seems 
uninteresting for several reasons. One is that it follows from (3.13) and (3.14) 
that it is impossible to match the shear stress ¢,,, which vanishes, with that 
given by (3.11) at some value of 7. Another is that it gives rise to a couple, 
associated with the shear stress tg, ”,, which is inappropriate for Poiseuille 
flow. 
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The Finite Flexure of an Aeolotropic Elastic Cuboid 


J. E. ADKINS & A. E. GREEN 


1. Introduction 


In a recent book by,GREEN & ADKINS (1960) a part * of Chapter II is devoted 
to the flexure problem for rectilinear and curved cuboids of aeolotropic materials. 
Some of this work is a special case of a number of general problems considered 
by ADKINS (1955) and, unfortunately, in re-writing the work for the book the 
wrong impression** has been given that some of the results for the particular 
flexure problem which is considered apply to a body possessing a completely 
general strain energy function W. If t’’ is the stress tensor referred to a curvi- 
linear system of coordinates #*, measured per unit area of the deformed body, 
and we take # to coincide with cylindrical polar coordinates (7, %, y;) in the 
deformed body, so that 

fee) Fah Va) (a) 


then t!2, 718 must, in general, vanish if the equations of equilibrium and boundary 
conditions are to be satisfied. In the problems considered by GREEN & ADKINS 
this implies some restrictions on the form of W (see below) and although the 
restrictions were in fact used throughout their work, some statements to the 
contrary are made, especially in §2.8. 

In a paper mainly concerned with general results in terms of a parametric 
representation, ADKINS (1955) examined the combined extension, shear and 
flexure of a cuboid and gave some results which, in principle, were valid for a 
general homogeneous aeolotropic body. We briefly re-examine this problem 
here in order to re-emphasize the conditions under which it may be regarded 
as solved for a material having a general strain energy function. 


2. Notation and formulae 
We quote the notation and formulae of GREEN & ADKINS (1960). The strain 
energy W of a homogeneous aeolotropic body can be expressed as a polynomial 


W = Wee;;) (2.1) 
in the components of strain e;;, where 
a oy" Oy" 


SISOS, DO, Psy Aiea. 
*xx We are indebted to Professor TRUESDELL for drawing our attention to this. 
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In (2.2) x‘, y' are the coordinates of points in the undeformed and deformed 
bodies respectively, referred to the same rectangular Cartesian system of axes. 
The stress tensor 7’’, defined in §1, is given by 


pe LE a a0 3 > 
a 2yI5 es Cea | OME ON oa ( 3) 
for a compressible body, where 
Ip=|2¢,, + O,6|- (2.4) 
If the body is incompressible, /;=1 and 
ij. 1 (OW aus oo o9) ij 
i 2 ee Ceey Or 07" nes (2.5) 


where # is a scalar function of # and G'? is the contravariant metric tensor 
of the curvilinear coordinates # in the deformed body. The equations of equi- 
librium, in the absence of body forces, are 


FiO, (2.6) 


the double line denoting differentiation with respect to curvilinear coordinates 
# in the deformed body. 


3. Flexure, extension and shear of cuboid 
We consider an undeformed cuboid bounded by the planes 


CA, Oe te =. Pe ee ew Si eh (3.1) 
We take # to be given by (1.1) where 
y—=71ces?, Ne.—7s5n9 (3.2) 


and we consider the deformation 
r=f(m), P=F2+y(m), y=Axytw(n), (3.3) 


where f, y and w are functions of x, and A, k are constants. Any plane x,=constant 
then becomes a cylindrical surface y= constant and the cuboid becomes a cylinder 
bounded by surfaces r=7,, 7, (4 >72) corresponding respectively to x;=A,, Ag. 
The constant A represents an extension ratio parallel to the x-axis and (xj), 
w(x,) represents shearing displacements. If m and w are non-zero the planes 
X%,= +B, x,=-+C become curved surfaces in the deformed body. 


If the body is incompressible then 
f=(2kx,+ Ky)! (3.4) 
where & is a constant, but y and w are unrestricted. 
With the help of (3.3) the strain components (2.2) are found to be 
26> 2 4 f* gles wil ein Ake 
2ey9 = f?/(A2 k?) — 1, 2e,=Aw', (3.5) 


2¢3=A7— 4, €y3=0, 
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and, from (2.4), 
I, = f? f'?/k?, (3.6) 


where primes denote differentiation with respect to x,. 
The stress tensor (2.3) for a compressible body has components 


fir anh ow 


722 — 


4 ow ! ko? ow yp’ ( a 
WAR SOGR: he) OG Afi Oe Come 


33. A2k OW kw’? ow Akw’ {( ow weeW, 
T , é a 7 a {= 7 ( T ] ’ 
ah 33 ih ery fel 0e3 6e34 
1 ow ow kg’ ow 
High ( 
2Af \ Ge. Olt fe,” 67) 
ices Ak (3 | new OW. 
2f \0as 6e31 if Ge,” 
Se ae 1 ( OW —-0W ky’ w’ ow 
2ff’ \ ess es» lain lara 
4 ki Pp (= ! ow w’ J ( ow e OW 
2ff \0as 0e31 2Aff' \ bey. . 0es, 


If the body is incompressible we must add terms #, p/r?, p to 11}, 2, 733 re- 
spectively in (3.7), and we can make use of (3.4). 

In the compressible case the stresses t’’ are functions of x, (or 7) only. The 
equations of equilibrium (2.6), in cylindrical polar coordinates (7, #, 3), reduce to 


dit Fide ap ae 
Che © Y : 
eee 312 
ay aA =0, (3.8) 
Hee “ile 
i cee 
ay ty o 


From the second and third of (3.8) we obtain 
Pe Ia eh, (3.9) 


where « and f are constants. If the surfaces y=7,, y=7, of the deformed cuboid 
are to be free from applied shear stresses, t!? and t%1 vanish on both surfaces. 


In this case 
p2=0, 13 —0, 


and « and f are zero in (3.9). Also, from (3.7), we then have 


, OW ow ow 
| = 
dS Ge, ' Oe, | Gea 3 


2w’ ow ow ow 


A Oey, 0e15 Gear 


(3.10) 
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If, in addition, there are no terms of the first power in é,, or ¢;3; in W then* 


one solution of (3.10) is 
gy =w=0. (3.11) 


We then recover the deformation considered by GREEN & ADKINS (1960) and 
all the formulae they give are correct ** if this restriction is imposed on W. 

Suppose W is not subject to this restriction. Equations (3.10), are then 
polynomial equations for the functions gy’ and w’ in terms of x, and 7(x,) and 
the first of equations (3.8) gives a differential equation for f(x,). When W is 
a general polynomial in the components of strain it appears to be very difficult 
to determine whether or not equations (3.10) yield an appropriate real solution 
for gm’ and w’; there may, of course, be more than one possible set of values 
of g’ and w’ satisfying (3.10). We assume, for the present, that at least one 
suitable set of real values of gy’ and w’ in terms of x, and f(x,) can be found to 
satisfy (3.10), and thus make 112, 71% vanish identically. The remaining stresses 
in (3.7) then reduce to 

eek Clie 
i ey ee 
1 oW ko? aw 
WRT Oey, ff > Berg’ 
33. 2k OW kw? aw 
ff’ @es5 fh dey” 

1 (ae aW\ ko’w’ ow 
2th \ Gees za fh e, ° 

Since the components of strain e;; in (3.5) are independent of x, and x; it 
follows that W may be regarded as a function of x,. Then, using (3.5) and (3.10), 
we have 


dw ow yer heey ees. ree 7] 
Te Se +H ot + PG 9” + w") + 


722 — 


(3.12) 


723 — 


Che eles ae ia ee pe \+ Aw” gs ow 
Obxy Ath? -* BN Oen) eeee Ak OUR 2 \ 05 ae 
ow ‘get , ow ite 
= an 2 Ran = ee 
Oba Ve te Cbs a wAURA 


and with the help of (3.12) and the first of equations (3.8) this further reduces to 


adW(%) 1 d 211 
Gh ok rag a Ms 


This equation can be integrated to give 


m= # [W(q) +H] (3-13) 
where Wy is a constant. From (3.13) and the first of equations (3.12) we have 

1p OW 

he ov =W(x%)+WM, (3.14) 


oi Terms containing the product e¢,3 are mot, however, excluded from W. 
This applies to §§ 2.8, 2.14. Corresponding restrictions on W are needed in 
the remaining sections 2.9, 2.10, 2.15—2.18 of Chapter IT. 
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a differential equation for /(x,), assuming that ’, w’ have been found in terms 
of x, and /(x,) from (3.10). 

When the body is incompressible the function /(x,) is defined by (3.4). The 
functions gy’, w’ are then determined from (3.10) and, assuming this to be possible, 
equations (3.13) and (3.14) are then replaced by 


= ee th = Win) +N, (3.15) 


a Loe a: 
PE Berg Seg bays? 
ow ow 
Saye "2 
Tt A a ae, ?, (3.16) 
che ale wl 7 di 1p BW 
Dh Ob, Cen 44° 


If the curved surfaces r=7,, 7, of the deformed cuboid are to be entirely 
free from applied stress then, in addition to (3.10), we must have 


W(A,) = W(A,) =—W, (3.17) 


for both compressible and incompressible bodies. Once again, further consideration 
is necessary before we know whether or not it is possible to satisfy (3.17). 


The faces x,=-+8B of the cuboid become curved surfaces in the deformed 
state. In order to find the resultant forces acting on these surfaces it is, perhaps, 
simplest to consider stresses acting across any plane #=constant of the deformed 
body. These are equivalent to a normal force N,, a shear force S, parallel to 
the y-axis, and a couple M, about the y3-axis where 


N= Jf 72722 dr dyg; 

Sa=JS[rrdrdy,, 

Me=/f[rrdrdy,, 
the integration being over the section #=constant of the deformed cuboid. The 
integrands can be regarded as functions of v and the integrations with respect to 


vy, can be performed immediately between the limits — AC + w(x.) and AC + w(%). 
Then, changing the remaining variable to x,, we have 


A, 
d 
N,=246 [ 7) am 
(3.18) 
=2AC (fees 
where we have used the first of equations (3.8) in the form 


f Peta (fat). (3.19) 
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A 
Also  F ye 
M,=24C [ fa dx, 

A, ; (3.20) 

=2aC {pays —k fT (m) + Wl anf. 
As 
If equations (3.17) can be satisfied then (3.18) and (3.20) reduce to 
NZ = 0; 

. (3.21) 


A, 
M,=2AkC ie Wahine AAW, 


These formulae, together with (3.18) and (3.20), are valid for compressible 
and incompressible bodies and are similar in form to those given by GREEN & 
Apkins*. They may be valid for a completely general strain energy function W 
provided equations (3.10) can be satisfied and, in the case of (3.20), provided 
also equation (3.17) holds. The doubt about their validity rests on our inability 
to decide whether or not equations (3.10) have suitable solutions for 9’, w’. 
By introducing the condition that W contains no terms of the first power in 
49 OF €3 (apart from é,, é,;) but is otherwise unrestricted, we have seen that 
one solution is y’ =w’ =0, and this is the case considered by GREEN & ADKINS, 
although they made some statements which gave a contrary impression. 

Forces are required over the curved faces of the cuboid which correspond 
to the plane faces y,=-—+C of the undeformed body, but we omit detailed ex- 
amination here. 

GREEN & ADKINS also considered the flexure of a curvilinearly aeolotropic 
cuboid and their work requires a similar restriction to be imposed on W**. 


* Chapter II, § 2.8 and 2.14. 
** In their notation W must contain no terms jointly of the first power in y49), Ya3)- 
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Propagation of a Plane Wave in an Isotropic Elastic 
| Material Subjected to Pure Homogeneous Deformation 


M. HAYES & R. S. RIVLIN 


1. Introduction 

In the present paper, we apply the theory [7] of the superposition of infini- 
tesimal deformations on finite deformations in an isotropic elastic material to 
the study of the propagation of a plane wave of small amplitude in an infinite 
body of the material which is subjected to a static, pure homogeneous deformation. 
It is seen that the secular equation for the determination of the square of the 
velocity of propagation in a given direction has three real eigen-values and 
correspondingly three real mutually perpendicular eigen-directions. Provided 
these three eigen-values are all positive, travelling waves may be propagated 
in the body in the direction considered with linear polarisations along each of 
these eigen-directions. If one or more of the eigen-values is negative for any 
direction of propagation, the body is inherently unstable in the state of pure 
homogeneous deformation considered. 

The condition for this to occur thus imposes a limitation on the strain-energy 
function which can be valid for a material that can be maintained in the state 
of pure, homogeneous deformation considered. In the case of unstressed materials 
obeying the generalized Hooke’s law, this has long since been pointed out by 
KELVIN [2]. It is quite apparent that a similar procedure can be applied in 
the case of visco-elastic materials, whether linear or non-linear, to obtain re- 
strictions on the constitutive equations. 

Finally, it is shown that the necessary conditions for stability of the material 
which have been derived imply certain inequalities proposed by BAKER & ERICK- 
SEN [3], although the converse is not necessarily the case. 


2. The fundamental equations 


We consider a body of isotropic perfectly-elastic material to be subjected 
to a finite static deformation, without body forces being applied, in which a 
particle initially at X,; in the rectangular Cartesian coordinate system x moves 
to x; in the same system. We shall call this state of deformation of the body 
the state A. 

The strain-energy function W, per unit volume measured in the undeformed 
state of the body, is a function of three invariants of the deformation J,, I, 
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and J, defined by a 
In = Gri 8s5 — 817 G3) (2.1) 
and 
I, = [es;|, 
where g;; is defined by 
_ 0% 0%; 
bt aan 2a) 


The components of stress o;; in the coordinate system x are given [4] by 


055 = sy [OG + LW) 8:5 — Bin Be We + LM 444], (2.3) 
3 


where W,, W, and W, denote dW/éI,, @W/éI, and GW/eI, respectively. Since 
the body is in equilibrium and no body forces are applied, 


003; 


=. (2.4) 


OX; 


We now consider that a further deformation, which may vary with time, 
takes place, in which the particle x; moves to x; where 


%;,= %;, + €U;. (2.5) 


é is a constant, which is small in the sense that terms of second and higher 
degrees in ¢ may be neglected in comparison with those of first degree. In the 
resulting state of deformation the body will be said to be in state B. Let 
o;;+€6;; be the components of stress in the coordinate system x, when the 
body is in state B, at the particle which is at X; in the undeformed state. It 
has been shown in the previous paper [7] that 


6,; = 215% {- 47,15 [¢;,(M41W,) — Bin &njg We — 13 We 6; 5] + 
+e (ht GM) Peles ew) = (2.6) 
aa RinOrg Vo Wa (8in Sj + Sin Se) + I3W, 6;;+, 


where 
= Ox; Ou; Ox; OU; 
Siig = aye 21.4 ~, (2.7) 
OX, OX, | Oks Oe 
1 =e Baz. 
a= Cieir— Cpee (2.8) 
= OU; Ox; OX, : Ou; Ox; Ox, . OU, OX; OXp \ 


1 
lee 21; e;jn 


(OX, OXy OX; "ON, Ok, Oy Wee eee ) 
and W,, Wy, W, are defined by 


WeEWi, i) (2.9) 
where se 


The equation of motion for the body in state B is 


a) a Aan 
Ox; (o;; + €%;;) =2@ one 
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' where @ is the density of the material in state B. If o is the density of the material 
in state A, then with the neglect of terms of higher degree than the first in e, 
we have 


@=o(1— ese). (2.12) 


Introducing (2.12) and (2.5) into (2.11) and neglecting terms of higher degree 
than the first in ¢, we obtain 


055,53 — F171 Un, 7 =O PUu,/0l. (2.13) 


The notation; is used to denote 0/6x;. In the case when the state of deformation 
A is a pure homogeneous deformation, o;;,=0 and equation (2.13) becomes 


F55,,;=0 Pur. (2.14) 


3. Small deformation superposed on pure homogeneous deformation 


We now consider the case when the state of deformation A is a pure homo- 
geneous deformation with extension ratios A, («=1, 2, 3) and principal directions 
parallel to the axes of the coordinate system x. The deformation is described by* 


A ekegs (3.1) 


Introducing (3.1) into (2.2), we have 


Bap = Aa Ay Sup = Aa Ons (3.2) 
With (2.2), we obtain 
3 
I, =F 2, We 
a=1 
3 
T= D dads, (3.3) 
“eB 


2492 42 
=A 243. 


Introducing (3.1) into (2.7), we obtain 


Ou Ou 
Bap = da ay + As ay 
a AS Cae SS A an 
= Ai Up y +Ag Ue, B- 
Introducing (3.4) and (3.1) into (2.8) 
eet 
I, ad > 21%, Uy a? 
a=1 
i 3 
To 2a a Aa) he (3.5) 


x The summation convention will not apply when Greek subscripts are used. 
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Introducing (3.2), (3.4) and (3.5) into (2.6), we obtain 


C= 2154 (W, + LW.) [— Un, I bap + NU, « so 1G Ux, p a 
3 
+ W, | 2 bap(23 Uy yt 2 ay u,,,) (A a, ae de »| + G6) 


4 Wy Iq tp, bap + [A2(0, + 1, W, — 2W,) + 15%] dah 


Thus 
Gap = 215 {M+ (Lh — 3 — 4h) Wa} A Me, + Ap Map), +B, 
and 


Bau = 25H, + LW) HE 2tt, 4 — a,8) + 
3 
+ We 13[ 22 0,» — Aaya) +25, By | + 
y= 
+ We Es Uy, + 
+ (+17, — 2%) + 147%) 


(3.7) 


We note, from (2.9), that 

AB (W, + I, W, — 2 W,) + TW, = {25 (Wa + ( — 42) Wor) + Lg War + 
+ {4alWie + (Li — Ag) Woe] +IsWe2} 12+ (3.8) 
+ {a2 [Mis + (i= 7) Wy] + I;Wys} Is. 


Introducing (3.5) into (3.8) and the relation so obtained into the second of 
equations (3.7), we see that 6,, may be written in the form 


3 
Gis, = 2 Gate: (3.9) 
p=1 


where c,g 1s given by 


Cup = 215 {[W, + Wa (Ly — Az — 45) AG (268 — 1) + AL ABW, + IW + 
+ 2 {Aa [Waa + (ly — 23) Wei] + Ig Wei} AB + 

+ 24a [Met (la — Aa) Woe] + 1g Wao} AB (Ly — AB) + 

+ 2(A[Ms + (Ly — Ad) Was] + Ip Wes} Js}. 


| 


(3.10) 


4. Propagation of plane waves 


We now assume that the displacements ew, arise from the propagation through 
the deformed material of a plane wave, with wave number , in a direction 
with direction-cosines /; in the coordinate system x and that w, is given by 


Ww; —— U; (t) a. (4.1) 


where 
a=, Xj, L1.= qe (4.2) 


Then, z is the distance measured along the direction with direction-cosines Lj. 


Plane waves in elastic materials 19 


By introducing (4.1) and (4.2) into the first of equations (3.7), we obtain 
Gp = 215 4*{W, + (I, — AZ — 23) We} (AR 1, Us + 131,U,) okett® (a +B). (4.3) 
Similarly, by introducing (4.1) and (4.2) into equation (3.9), we obtain 


Coa ace Bo tne (4.4) 
y= 
Introducing (4.3) and (4.4) into (2.14), we obtain an equation of the form 
EW opt Mo de 07/82?) Cais (4.5) 
where Sap = Ag lala (A - 26,9) [Wi + We (I, — 42 — AR) + (4.6) 
+ Sap] DAE + Walls — 2 — )]} + EB cap lal. : 
Writing 


Cag = Cag + 215 4{[W, + We (Ly — Ae — Ap)] Ag(25ap— 1) + AGARWe + IgM}, (4.7) 
we obtain from (4.6) 


Sap a (Ag 3 W a I, W3) lls a 


3 (4.8) 
+ bxp{ LBB + Wa (h— 22 — 49]} + $1 Cap lap. 
y=1 
Replacing 0/0t by tm in (4.5), the secular equation may be written as 
‘ Isep— C*Sag] =O, 49) 
ae C2= 11h 9 w2/h?. (4.10) 


It follows from (3.10) and (4.7) that ¢,, is symmetric. Hence from (4.8) s,, is 
also symmetric. 
} 5. Discussion of the results 
Equations (4.5) and (4.9) can be transformed to any other rectangular 
Cartesian coordinate system y, related to the system x by the transformation 
Vi = 455 Xj, (5.1) 
with 
Wp jp = Ani 4,5 =9;;, [a,j] =1- (5.2) 
Let <V,(¢) e'** be the displacement components in the system y. Then 


We note that z is unaltered by the transformation from the system % to the 

system y. With the relation (5.3), it is evident that equation (4.5) can be ex- 

pressed in the form (t;, + Iho 8, 2/0) V;=0, (5.4) 

where 

bij = jp Aq Spq- (5.5) 

With the same notation, the secular equation (4.9) can be written in the form 
2* 
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Since ||s,g|| is symmetric, the coordinate system y can be chosen in such a way 
that ¢;,,=0 (¢=-7). Then, it is seen that equation (5.4) takes the form 


173 Q2 0 
(teat gat Gp) a= a 


Thus, the directions of the axes of y are the three directions of polarisation for 
which waves, which are not mutually coupled, may be propagated in the direction 
of z. The velocity of propagation, along the direction z, of the wave polarised 
in the direction of y, is (t,4/31%0)?. t,. («=1, 2, 3), are, of course, the solutions 
for C2 of the secular equation (4.9). Provided all of these solutions are positive, 
the velocities of propagation of all three waves are real. 

However, if one or more of the solutions of (4.9) for C* is negative for any 
direction of propagation of the plane wave, 7.e. for any real values of /; satisfying 
the relation /;/;=1, it is easily shown that the material is inherently unstable 
in the state of pure homogeneous deformation considered. Equation (4.5) may 
still be written in the form (5.7). Let us suppose that 4, is negative, then the 
general solution of (5.7) with «=1 is 


GSA ase, (5.8) 
where A, and B, are integration constants and 
Pe as 
Oo, = k (— beat EE P “ 
: 313 oe 


The displacement component « J, (f) e‘** is then given by 
eV, (Dew t= eApeey Re Byer (5.10) 


We shall now assume that when ¢=0 the rate of change of the displacement is 
zero. Then, from (5.10), 4;=B, and we can re-write (5.10) as 


eV, (t) fF Remco aS yeaa let’ (5.11) 


It is apparent that a displacement parallel to the axis y,, which is sinusoidal 
along the direction z with wave number k, builds up exponentially with time. 

We thus have the result that if one or more of the values of C? given by 
(4.9) is negative, for any direction of propagation of the plane wave, the material 
considered is inherently unstable in the state of pure homogeneous deformation 
described by (3.1). 

For a real material which can be maintained in the state of pure homogeneous 
deformation (3.1), the dependence of W on J,, J, and J, must therefore be such 
that the three roots of (4.9) for C? are all positive, for all J; in the expression 
(4.6) which satisfy 1;/;=41. The conditions for this to be the case may be written 
in a variety of forms. For example, a necessary and sufficient condition that 
the roots of (4.9) be positive is that the quadratic form S;;¥;¥,; be positive definite. 
Alternatively, the necessary and sufficient conditions may be written 


trS>0, (XtoS}*-itrS4>Omand |She.0; (5.12) 


where the notation S=||s,,|| is used. 
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We note that the conditions (5.12) are independent of frequency. Consequently, 
if they are violated for some /; and for some frequency, or wave-length, they will be 
violated for all frequencies and-wave-lengths. If the conditions are not satisfied 
for some values of J,, J, and J; and for some real value of /;, satisfying /;1;=1, 
any infinitesimal vibration present, due to Brownian motion or other cause, 
will build up. The material cannot then be held in the state of pure, homogeneous 
deformation considered. Conversely, if a material can be held in the state of 
pure, homogeneous deformation considered, then a strain-energy function which 
violates the condition (5.12) for this state of deformation and some value of / 
is inappropriate to the description of the properties of the material. 


6. Relation to the Stability Criteria of Baker and Ericksen 


BAKER & ERICKSEN [3] proposed as a necessary restriction on the strain- 
energy function for an isotropic elastic material, that, when subjected to a pure 
homogeneous deformation with principal extension ratios A,, the relation 


W, + 2,W, 2 0 (6.1) 


must be satisfied. It is easily seen that the relations (5.12) imply the relation 
(6.1), while the converse is not the case. 


We may re-write (4.8) as 


3 
Sup = Oup Qu Aay li + (1 — b42) Byplals, (6.2) 
ay | 
where 
pee han Ove Marea L147, 5) + AB{W, + (Lh — 43 — 23) Wo} ae 
an + : 
Bup = 19 Ag We + IgWa t+ $18 Cup (« = B). 
Consequently, 
3 3 o 2 3 > 3 
Up pes ad te ool ( Y Aaals) = 2 » Aus] 
c= e=L p= at a=1 
and 1D fu (6.4) 
trs?= > | SA ea) Beh El 
yy. Se) 
From (6.4), we obtain 
3 
(tr S)?— tr S?= intr hns (6.5) 
Ys n= 
where 
3 
U,, = D> AeA Aa Ae) ao (1 =S057) Bee (6.6) 
“a, p=T 
The conditions (5.12) may therefore be written as 
ae 
p=1 tA 
3 
DU nly = 0 (6.7) 
yn=l1 


and 
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These inequalities must be valid for all values of /?, /3, 13 which satisfy the relation 


ee (6.8) 


a=1 


Necessary and sufficient conditions that the first of the relations (6.7) be valid are 
age: (6.9) 
1 


A necessary condition that the second of the relations (6.7) be valid is that 
Uz > 0. (6.10) 


With (6.6), the relation (6.10) may be written as 


DAA, eso. (6.11) 

B,y=1 

By 
Finally, a necessary condition that the third of the relations (6.7) be valid is 
Ay, Ag, A3q > 0. (6.12) 


From (6.9), (6.11) and (6.12), it can easily be deduced that 
Ag 20: (6.13) 


The inequality (6.1) derived by BAKER & ERICKSEN is a special case of (6.13). 
For, from (6.13) and (6.3), we have, if «+6, 


WA 2 A) 0. (6.14) 


Bearing in mind the definition of J,, it is evident that (6.14) implies (6.1). 


It is apparent that, in general, the relations (6.1) do not imply the relations 
(5.12), since the latter relations involve the second derivatives of W with respect 
to the strain invariants J,, IJ,, J;, while the former do not. 
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On the Motion of Two Linked Bodies 


C. MASAITIS 


Communicated by E. LEIMANIS 


1. Introduction 
The use of satellites for weather observation, communication, navigation, 
| etc., requires a certain orientation of these bodies in space. Various devices 
/have been proposed for the purpose of damping the nutation of a force-free 
body and thus stabilizing its orientation. Two examples of such devices are 
discussed in [7] and [2]. However, the complexity of the equations and the 
limitations on the range of the variables make the complete analysis of their 
motion difficult. Furthermore, the shifting of the center of mass within these 
systems has obvious engineering disadvantages. 

A system of two linked bodies, proposed by S. ZARoopNY, seems to offer 
‘the same possibilities as the two systems mentioned above but is free of their 
disadvantages. In [3] he and J. W. BRADLEY have exhibited some numerical 
solutions which show that with very large friction of a viscous type, 7.e., pro- 
portional to the relative speed of rotation of two bodies, the motion is of the 
form of a perturbation on the free nutation. Moderate values of the coefficient 
of friction produce a damping of nutation or metastable forms of nutation, 
depending on the initial conditions of motion. Very small friction produces 
rather complicated motion with markedly different position of points of equi- 
librium from the case of rigid body. In some cases of the motion in the vicinity 
of the points of equilibrium it was difficult to distinguish between the accu- 
mulation of numerical errors and the metastable character of the motion. There- 
fore, it was desirable to attempt the analytic treatment of the problem. 

Although no elementary solution can be constructed, the equations of motion 
are sufficiently simple to offer possibilities for a perturbation treatment. The 
unperturbed system can be chosen as the one without friction, which admits 
a solution in closed form. In this paper such a solution is constructed and the 
nature of the equilibrium points investigated. The results could be used for 
selecting special configurations of a case with friction, for investigation by 
numerical solution. 


2. The Equations of Motion and the First Three Integrals 
The system to be considered consists in two rigid bodies B, and B, mounted 
on a weightless rigid axle A passing through the centers of mass of B, and By. 
Let O be the center of mass of the system, J,, I,, J, be the principal moments 
of inertia of B, about O, not all of them equal, say 


(1) dae 
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Let the principal moments of inertia of B, about O be J, J, and J,. Suppose 
that the axes of the moments J, and J, coincide with the axle A, and that B, 
and B, can rotate about A without friction. Let the coordinate axes (1, 2, 3) 
(Fig. 1) through O be parallel to the axes of the moments J,, J, and I,, respec- 
tively. Suppose that no external forces or moments act on the system. Then 
QO can be chosen as the origin of an inertial coordinate system (x, y, 2). If 
@ = (@, M2, 3) is the angular velocity of B,, where w,, 2,3 are components 
along the axes (1, 2, 3), M the moment of force acting on B,, and H,=(I, ,, 
3 I, @., I3@3) the angular momentum of By, 

about 0, then, since O is the origin of an 

inertial coordinate system, it follows that 

HN where 4 = 44x anda 

is the derivative relative to inertial axes 


| (a, V2), = is that relative to rotating axes 


(1, 2,3), and @ is the angular velocity of the 
latter. Therefore, the equations of motion of 
pote I, @,+ (3 — I,) 203 =M,, 
(2) I30@+ (1, — Is) 30, = Mb, 

I3@3 + (12 — 1,) @, @, = M3. 


Because of the assumed axial symmetry 

of B, the angular momentum of B, about 

O is Hy=(I 1, I wy, 1,3), the components 

being taken along the axes (1, 2, 3). Let the 

Fig. 1 moment of force about O acting on B, be 

M’' = (My, Mg, Mg). Since B, is constrained to 

rotate on the axle A, there follows w;=@,, @3=@,. Let w;=«w,. Applying the 
operator d/dt defined above to H, yields the equations of motion of Bag 


A 


To, + I, W,@,4 = I W3 Ws = M,, 
(3) I @, — 1,@,0,+ 130, = M3, 
TiO me Be 
Since no external forces or moments are applied to the system, one has 
M=—M,,, M,=—M;. Furthermore, since both bodies can rotate independent- 
ly about axis 3 without interaction, it follows that Mj=Ms=0. Hence (2 ) and 
(3) yield 
(+ 1) @1 = 2 [(I, — Ig + I) wg — I,@,4], 
(Lo = Lo oy ee a I) ws ; Te@ i, 
I3@3 = @2(I, — I), 


TO, ==0% 


(4) 


Clearly I,@,4=C, a constant. Let the direction of the axis 3 be so chosen that 
CeO: 
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Let Pas ieee! 
nes STI 
Se daietly 
a= Beh ae 
b= 757 >0, 
b= oz > 0. 


By substitution of this in (4) one gets 
Oy = @2 (4,03 — by), 
(5) Wz = Wy (az 3 + by), 
Ws oe a ag Oy Wo . 
Multiply the equations (5) by a3 b.@,, d3b,@., (a, 03+ a, b,) w3, respectively, 
add, and integrate. This yields 
(6) az by @} + ag b; 03 + Bos =B f?, 
where 


B= a,b, + a,b, >0, 


and #? is a constant of integration. Similarly, multiplication of (5) by a, a4,@,, 
— A, @,,, —f, respectively, followed by addition and integration yields 


(7) Ay Ax Ot — A, A,03 — 28 WO, = 28 9, 
where q is a constant of integration. 
3. Singular Points 


In the examination of singular points, 7.e., points of equilibrium of (5), it 
is convenient to consider the following cases: 


(a) T=i,—I,, 1.é., a<—0,— @5 > 0; 
(b) l=I,—T,, ICs eae O Ce ae OF 
(c) Tel eed dd 1.C2 a 0 Ag 05 
(d) Pt, 5, 1€., a >0, a=0; 
(e) I>I1,—1,, ies Vd, eh dan 
In cases (a), (c), and (e) the sets of singular points of (5) are defined by 
(i) oath aes Bs 
b. 
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In case (b) the set of singular points is defined by (i) and (iii) above, in case 
(d) by (ii) and (iii). It is of interest to examine which singular points in the 
phase space (@,,@y ,@ 3) are points of stable equilibrium and which are not. 
Let U(P, «) denote a spherical neighborhood of the point P with radius «. It 
is said that the point of equilibrium* P is stable if for every ¢>0 there exists 
a 6(é)>0 such that if the motion starts at any point in U(P, 6), it will remain 
in U(P, e) for all values of the independent variable ¢. 


The solution of (6) and (7) for a,j and a,5 in terms of a, yields 


(8) az, = 4, (p? — ws) + 20,9 + 2b, 
(9) Az 03 = A, (p® — w3) —'2b.q — 2bg0s. 


Lemma. Suppose that p and q are so chosen that (8) and (9) are satisfied by 
the coordinates of a point of equilibrium P. If with these values of p and q (8) and 
(9) define a non-degenerate curve passing through P, then P is a point of unstable 
equilibrium. 

The possible types of equilibrium points of (5) are: (1) nodal points, (2) focal 
points, (3) centers, (4) saddle points. However, any solution passing through a 
singular point P is an intersection of S, and S,, the quadric surfaces defined 
by (8) and (9), respectively. The intersection has at most four branches passing 
through P. Hence P cannot be a nodal or a focal point. If P(@,, @,, @3) were 
a center, then w,=@,, ®;=@ , ®3=; would be a solution of (5) consisting of 
a single point, 7.e., a degenerate curve. Hence the hypothesis of the lemma 
implies that P is a saddle point, and hence it is unstable, g.e.d. 


Case (a). Points of unstable equilibrium. Let A be an arbitrary point on the 
line (ii), say, A(o, By, i). In what follows the constants of integration cor- 
responding to initial conditions (@,, @2, 3), (@,, @,, @3), efc., will be denoted 
by (f,49), (6,9), etc., respectively. Since only the second power of f appears 
in (6), (8), and (9) it will be assumed that sees By substituting on (0, é,, 2) 


ay 


in (6) and (7) it follows that ¢ = — a (8 4). With this value of ¢ formulas 
(8) and (9) can be written as 


2 
(10) a, 004 = — a («0,— 21) , 
ay 
(11) A, (03 = — Ay 3 — 20, 0g B p> iO 
by ay 


If D is one pele of the discriminant of the right-hand side of (11), then 
a 
alt : (224 | nt). Since in this case a,<0 and B=a, b, +a, b,>0, it follows 


that ~hs8 and hence au Rlie Dy ba 


ay As Ay As 


b, b, 
pe , hence p+ ray a >0,1e., D>0. Therefore w, as defined by (41) is real, 
L2 


Since ess (6) implies that 
1 


* Stability in Lyapunov’s sense. 


: Let fh <t,<— a 


Motion of two linked bodies 27, 


| provided that 


(12) pain le eee, 
Ay a a Ay 
| bob enc eee ae pe 
N F=2(4 12); A Peal ieee 
Ow Be a = a therefore a <(p?+ a a or a3 <D, whence 


| —D< + <yD, Le., —|D<b,+ 2" <)D. Therefore — A aD < a < re SBN 
1 1 


ay a, 


j 2 
| Since w, as defined by (10) is real for all values of w,, this last inequality together 
} with (12) shows that A lies on a non-degenerate integral curve of (5). Hence 


by the lemma, A is a point of unstable equilibrium, 7.e., the line (ii) consists of 
points of unstable equilibrium. 


Next let B be a point of the line (iii), say, B(0, 0, @3) (®;<0). The constants 


of integration are s=G=—@,. With these constants (8) and (9) reduce to 
| (13) az 01 = (3 — Bs) [(20, — a, ®) — a, 0g], 
(14) dz 03 = (3 — Oz) [(— 2b, — ay Wy) — ag og]. 
by 


= Then fy< — 222 — Gy. Let I be a closed interval |@z, 
2 2 
As ~ 


| ——#—@,|. Both factors of (14) are non-negative if w,€J, 1.e., ws is real for 


As 


2b, 


all values of 3 IN a Again since i <O; it follows that =1—8< 3. Hence 


1 
a or 2b, — — A, 3 > 4,3. Therefore for every value 


of ws in J, w, is eal Since @3;€J/, it follows by the lemma that B is a point of 
unstable equilibrium. 


Points of stable equilibrium. It follows from (8) and (9) that for e,>0 there 
exists an 7>0 such that 


(15) PP la, gama, 
and 

(16) | ws — @3|<n7< & 
imply that 

(17) Vee aon ee eee 


In fact, by (6) and (7) for every y>0 there exists a 6>0 such that 


(18) |@; — @;| <0, CHAO 
implies 
(19) Bot ae | ae 


If (w,, 2, @3) is a point on an integral curve of (5), then w, defined by (9) must 
be real, and therefore 
(20) eae ee 


As a, 
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where R2—R2(p, q)=03-+ a3 p? — 2a, b.q; R,=0. Also, (20) can be written in 
the form 


(21) 


b 
3 + — 
OE 


: ; ; ~ b 
Let E be an arbitrary point on the line (i), say, E(,,0, — 74). By sub- 


2 
stituting the corresponding constants of integration #? and 7 in \k, it as\seew 


that R,=R, (p, 7) =0. Hence for every 7>0 there exists a w>0 such that 


(22) IP-Bl<n, |¥-al<n 
imply 
(23) Ro(f, 7) < 420- 


Let e>0 be arbitrary. Since for @,=0 E is a point of the type F, which is 
min( 5 a. ore € ie 
2) 


discussed below, it can be assumed that @,==0. Let e.= agi 


and choose 7 in (15) and (16) such that (17) is satisfied. With A so defined 
choose wz in (22) such that (23) is satisfied. In (19) let ymin (mw, 7), and choose 
6>0 in (18) such that 6<|@,| and that (19) is satisfied. Now let Q(@,, @,, @s) 
be an arbitrary but fixed point of U(E, 6). Then (18) is satisfied, and hence 
(19) holds with #2 and q being the constants of integration for the motion starting 
at Q and tracing a certain curve J”. Since (19) holds, it follows by definition 
of » that (22) is satisfied; hence (23) holds. For all points of /’, #2 and g have 
the same values, namely # and g. Therefore (23) and (21) imply (16); (45) is 
implied by (19). Thus, (17) is satisfied for every point of J. Since ®,=0, (17) 


ey 1 
on ie Or, 


implies that w2<—, or |w,| =|@,—@, <E. Also, by (17), 
2 9 2 2 y 


ie., —$@i<w{ — Wi, or 403<w?, ie., wo? +0 for all points of I”. Since |@, — 6, <a 
6<|@,|, ®,>0 implies that —@,<@,—@, or ©, >0; similarly ©,<0 implies 
that @, —@®,< —@, or @, <0, 1.e., @, has the same sign as @,. Since @, starts from 
@, and varies continuously on J’, and since w7=-0 for every point of I’, it follows 
that w, has the same sign as @, for every point of J. ps |o, + @,| >|o,| 
& 
; |a,+ @,| 
definition of ¢,, |w,—@,| <— on J’. This completes the proof that I’c U(E, e). 


for ever point of J. Hence by (47) |a,—S,|< << a Therefore by 


Hence the points on the line (i) with w,==0 are the points of stable equilibrium. 
Let f° be a point on the line (iii), say, F(0, 0, @). With the corresponding 
new interpretation of @,, ®,, ®s, p?, J and arbitrary e>0, let e¢,=min (<, ry 


and choose 7 in (15) and (16) such that (17) is satisfied. Now #2=@2, ¥=—a,, 


and R2= (b, -++ 42 @;)*. Therefore there exists «>0 such that if (22) is satisfied, 
then 


(23’) | Ro — |bg + ay @gl| < ayn. 
Let y=min(, 4) and choose 6>0 in (18) such that (19) is satisfied. Suppose 


now that the motion starts at some point (@,, @), ®,) € U(F, 6). It follows from 
(7) (since a,<0) that for every 420 we have w,;=—Q; hence 3 —®,>—G+q. 


: 
a 


| 
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By definition of 6 the last inequality can be written as 


: (24) 3 — Ws > — 7. 
Suppose now that 
(25) es, 
2 


e., b,+a,@,;=20. Then it follows from (23’) that = ae 


by (20), w;—@3<~y for all t. This, together with (24), sigs that (16) is satis- 
fied. Since (15) is implied by (19), (17) follows. Since now @,=@,=0, (17) 


implies that w?< a or |w,| =|o,—4,| <i. Thus, F is a point of stable 


2<@,+%7. Hence, 


equilibrium, 7.e., all points of (iii) for which w= = are points of stable equilibrium. 
2 


Next consider the case Oy a zi 1<0. Let Ri(p, g)=)3+ a? p24 2a, bg; then 


R? = R?(B, 9) =(b, —a,@)2>0, or ioe a, 3 —b,>0. There exists “> 0, such that 
for every # and J, satisfying (22), R?>0 and 


(26) |R, — R,| < |a,| 7. 
With this new definition of ~, and with 7 ie as above, again let y=min (w, 7). 


' Now choose 6>0 in (18) such that b< a, and such that (19) is satisfied. 
_ Then (22) and (24) hold; (22) implies that R?=R?(6,7)>0. [f and 7 are as 


before constants of integration corresponding to Q(@,, @,, @3) € U(F, 6).] There- 
fore it follows from (8) that for every point of the integral curve J’ passing seueh 


Q either oe sears 


Or Ws ria But certainly Dy<Oy+0< Gi, + 2) — Gy. 
Since — 7 503 this implies that Bye Aa Since @3 corresponds to s point 


of J’, this implies that ost, Since w; varies on [’ continuously, it follows 
that : 
(27) aS Es 

ay 


for every point on I’. Replacing (26) by Ree R,<|a,| 4 or 4,03 —0, — R,<|a,| , 
, and substituting from (27), we see that 


(28) Ws — O3< 7. 


Now (24) and (28) show that (16) is satisfied for every point of J’. Since (15) 
is implied by (22), (16) holds for every point of J’; thus it follows again in the 
same way as above that IC U(F, «), i.e., the points of (iii) for which w3< fai 

1 
are points of stable equilibrium. In summary of the results above it can be stated 


, ‘ 5 b 
that the points (or, 0, — +1) and (0,0,@,) with w3= — _ or Wg<——k are 
2 


2 1 


uniformly stable for every k>O. The points in a neighborhood of (0, 0, a) with 
ay 


b 
toy < Ft ave stable, but not uniformly. The points (0,0,@ 3) with ae 203 a 


and the points (0, Wo, 2) are unstable. 
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: : - nae b 
Case (b). Unstable points. If a point B is on the line (iii) and a3;< arr 
2 


it follows in the same way as in case (a) above that w, is real provided w3€J. 
Since now (13) reduces to a,@{=@3+ p, it is seen that @, is real if w3€J, thus 


the points (0, 0, @s3) with wz3< — 22 are unstable. 
2 


Stable points. The argument of Case (a) can be repeated to show that all 
points of the line (i) and those points of (ili) for which w= 2 are also stable 
equilibrium points. 

Case (c). Unstable points. Let B have coordinates (0, 0,@ 3). The solution 
of (5) is again given by (13) and (14). If @< —22, it ate: as before, that 


2 
for every @s; se T, @, is real. ara now os > 0, it follows that 2 > - Os 
ay 2 
—@3>@3 whence it follows that @, as defined by 


ay 
(13) is real. Therefore, by the lemma, B is an unstable point. If 3 > —— i i>0, 


Her; De 


then 2 = —@3<@,. Hence if c=max (0. 2 A ,), it follows that c< o, and 


1 af 
therefore a closed interval I’=[c, @,] exists. If w,€J’, then the first factors 
of (13) and (14) are never positive. The second factor of (14) is negative for 
all values of w,20. Hence w,€J’ implies that w, is real. Now, for w,€1’ we 


have O32 2 a —@3; hence the second factor of (13) is never positive, 7.€., @, 
is real yorued w3¢J1’. Therefore the point B is again unstable, by the lemma. 
Stable points. Consider E (a, 0, —2), i.e, E is a point on the line (i). 


2 
Since (20) holds also in this case, the argument about the line (i) in Case (a) can 


be repeated and the same result obtained. If G(0, Os, 2) is a point of the line 


lod ~ . . 1 . YY 
(ii) and #, 9 are the corresponding constants of integration, then R, (f, 7) =O. 
Furthermore (8) now implies 


(29) be aoe, 
ay ar na ay 


The argument of Case (a) about the line (i) can be repeated using this last 
inequality to lead to the conclusion that the line (ii) consists of stable points. 


Stable points on the line (iii). Consider a point F(0, 0, @3) with —?s <0, a Oy 


Since (20) holds in this case as well as in the Case (a), it follows as in Caseya (a) that 
(30) 3 — Og <7. 


Now 0, —a,@,=0; hence R=) — a, , and therefore it follows from (26) that 
- a ©; > —y, and therefore by (29) that w;—@,>—y. This, together with 
(30), establishes (17), and the argument can be completed in the same way as 
for the point F with ®,=> — 2 in case (a). Thus, in this case the points F(0, 0, ®,) 


x b 
20.5 fa are stable. 
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| The transformation w;= 3, @2=1, ®;=Ws reduces Cases (d) and (e) to the 
_ Cases (b) and (a), respectively. The complete classification of the singular points 
_ is given in the table below. 


Table. Equilibyrium Points 


Case Stable Unstable 
b, 1 
(a) =<, <KC9, On=O, OF ieee @,=0, —©O<a,<0, o,=— 
2 ay 
b b 
2 
@,;=0,=0, #,2—— C=O, = OF = a, 
Ay a 
b 
= 1 
W,;=@,=0, On 
a 
b bs (Oe 
(b) =C< MO Ke, Os—=—O, Oy= ®,;=W,=0, &3<—— 
ay Ga 
bs 
=O, Wy 
Ds 
by 2 
(c) SSK, O=O, Gy =— O,=W,=0, W3<—— 
Ay ae 
by j by 
O,=0, —W<@,<~, W3s=— WO;=W,=0, > — 
ay ay 
b b 
@1==W,=0; 2 <0, 5 + 
ay ay 
d by by 
(d) 0,=0, —oO<a,<0, o,=— O,=0,=0, > — 
ay ay 
b 
O,=W,=0, W3S—+ 
ay 
by be 
(e) ®,=0, —~©O<0,<0oo, aaa —0o<@,<00, W,=0, Oe 
1 2 
b b b 
O,=0,=0, o> —— @,=0,=0, —<o,<—— 
Ay ay a, 
by 
Op =O5= On OSS 
ay 


4. Complete Solution of the Equations of Motion 


Substitution of (8) and (9) in the third equation of (5) yields an equation 
in ws, which can be reduced to quadratures. The solution is expressible in terms 
of elliptic functions (including the degenerate cases with modulus x=1 or 0). 
The exact form of a solution depends on the character of the zeros of the right 
sides of (8) and (9), i.e., on the sign of R{ and R}. Since Rj and Rj are linear 
in g, it suffices to examine the signs of Rj and R} at extremal values of g for 
any fixed value of p. A simple calculation shows that the extreme values of q are: 
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b 
(31) “—h=5 f = (6% a a 
(32) am (6 
(33) a i tae 
(34) SF) 
(35) Q2 alte 
(36) Ri(q) = 45 ap (62+ 2 | 
(37) Ri(q2) = 0, 
(38) Rif qs) = (b, — a, p)*, 
(39) eT ee 
(40) RG) = 
(41) Ri(q2) = wf (oe ete), 
(42) Re (93) = (bg + 42 P)?, 
(43) es (4s) = (b, — a, p)*. 


The zeros of the right-hand side of (8) are 7, and 7,; those of (9) are 7; and 7,, 
where 
Psd Pa pe ts VRE = 0 Pees ripe Ho b,— /R3 


= to = 
) 2 ? 3 ? 
ay ay a a 


Case (a). It follows from (6) and (7) that a solution of (5) is real only if 


(44) —PS4. 
But (9) and (40) imply that @, is real if and only if 
(45) q<H- 


Since a,< 0, it follows from (8) and (9) that m, and w, are both real for properly 
chosen values of w, if and only if 


(46) 135% 
or 
(47) Ue) = V4 , 


provided 7, and 7, are real. (47) is satisfied only if 


(48) so and |g|<, 
or 
(49) q> 72 and p> 


2 As 


2 aes 
_ which is equivalent to abieb EVR ee AT 8 By adding a 1 p2 2bs4 
ae 
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The second inequality of (49) follows from the first and (45). If es then 
} 2 

| prt “1%2 <9 and hence by (31) ¢.<49,; by (37) and (45) Ri(q)=0, ie., 7, and 
1 5 hed’ — 
4 7, are real. By (45) = Since in this case a, g — b, <0, one has a, q¢—b,<\R? 


a : 
A, Aa, a, ag a A a, 


_ to both sides of the above inequality one gets 


2 b3+a3p?—2a,b,q , 2BVR3 _ b2+a2p?+2a, 6,9 
< ek 1 5 ihege | 5 
ajag * a T aa = a ; 


that is, 


ey REG Rs Pees 


a Ag a ay ay A Ay ay 


By subtracting 0,/a, from both sides of this inequality one gets —7,< —17. 


Therefore by (46) and (47) in view of (48) and (49) it follows that |g] <p. Hence 
by (39) and (43) Ri(g)>0 and R}(q)>0, i-e., 73-74, unless Ee When 


2 
—p<q<hf, all roots are distinct, and 7,<7,. Then w, is in the form 


_ Asn2(C,t+C,,*)+B 


60) A Cx8(C. CD” 


where C, is a constant of integration and A, B, C, D, and C, are known functions 


of 4%, 72,73, and 7,. The form of these functions depends on the position of the 


initial value of w, relative to 7,, 72, 73, and 74. If |g| =, then the point (0, 0, @,) 
with |@,| = lies on the curve defined by the corresponding solution of (5). 


Since pee or @, = — 22, by the preceding section (0, 0, @3) is a stable point, 
2 2 
and therefore by the lemma the solution is a degenerate curve consisting of a 


single point (0, 0, @s). 
Next consider the case 3 =p | —— Then by (31) ¢2=9,, and hence 
2 12 


_ by (37) Ri(q,)=0. Rj is a decreasing function of g. By (44) and (45) —p<q<q,; 


hence Rj (g,)=0 implies that Rj (7) =0, while (40) implies that R3(7)=0. There- 
fore 71, 7%, 7%, and 7, are real and distinct except when g=q,, 93, @4. Thus, when 
G+ 91> %3> Ia» Wy iS of the form (50). 

If g=43, then the point (0, 0,@3) with @,;=£ lies on the curve defined by 
the solution corresponding to the constants # and q,. Since ®=—q,=p, it 


follows that eee hence @3> — 22, and the solution is degenerate. 
a 


2 2 
If g=q, then the coordinates of a point on the line (i) satisfy the solution 
and again the solution is degenerate, as follows from the lemma, since the line 
(i) consists of stable points. If g=q,, then 7,=7,—= — , and 


(51) sz = + (3 + P) (20, +a, 2 — aa) ( 2b, + dap — Ay Qs). 


The integration yields 
__ Atanh?(C,7+C,)+B 
(52) 3 = Ctanh?(C,?+C,)+D ’ 


where C, is a constant of integration and C,, A, B, C, and D are known functions 
Dl 7) 4,;\and Dp. 
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If a _ bb x ah ay then by (31) ga<q%-. Since, as before, the range of g is 
Ay 


between — aa g,, it follows from (37) that Ri(q) changes sign when g=4q. 
By (40) R3(q)=0. Therefore for g<q, and q not one of gg, q4, all roots 741, 72, 73, 
and 7, are real and distinct. Accordingly, this leads to a solution of the form 
(50). If g=gg, it follows again in the same way as above that the solution is 
degenerate. The same holds for g=q. If g=q,, the argument leading to (52) 


su es b 
requires no modification, i.¢., 3 is given by (52). If g=q., then 4=7,=—~. 


ay 
Therefore | 
: b : 2 B b, b 
(53) @, = (0 a |/ ay ( Ay 03 — 25,03 - 5 p24 aah 
If p< = integration of (53) yields 
if 
A cane (Crt Cs) 
(54) Oa — Cte (Gite Can 
F See : b 
where A, B, C, D, C, and C, have meanings similar to those in (52). If p= — = 
then (53) can be written in the following form: : 
: b, \$ 
(55) 3 = + (05 — a |) a2(o5 + 2 28 mak 
Integration of (55) yields 
_ AC t+C,)2+B 
(50) O2=< OC Oa =E Ds 
If po, then the integration of (53) leads to (52). If qa<q<q, then all 


roots are distinct and R?<0, i.e., 7, and 7, are imaginary. In this case 


_ Acn(C,t+C,)+B 
(67) 8 C oniC, ECan 


Case (b). Here @§ is a third-degree polynomial in w;. (9) implies that the 
solution is real only if all three roots of this polynomial are real. The roots are 
distinct if g is none of 9,, 93, q4, and in this case the solution is of the form (50) 
(for certain regions of the initial pues of ws one of the constants C or D is 


ZeTO). If G=9,, Ge OF g==g, and ps2 , then by using the lemma and the results 
of Pipes section it can be ee that solution is degenerate. If g=q, and 
pe] Pa , then w3= A tanh?(C,t+C,) which is of the form (51). 


ee (c). Since in this case a,>0 and a,>0, then (8) and (9) show that for 
a real solution all four roots 4,72, 73,7, must be real. For fixed # now R? is 
an increasing function while R} is a decreasing function of g. Hence it follows 
from (37) and (40) that ¢,<q<q and that 7,=7, only if g=q, 44=7, if g=q,. 
A simple computation shows that any other combination of equal roots ee, 
only if g=q, or gy. Therefore if g+q;, i=1, 2,3, 4, then Ws 1s given by (50). 
Y= OF gz, then by using the lemma and the fact that the lines (i (i) and (ii) ee 
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‘of stable points it can be shown in the same way as above that the solution is 


| degenerate. The same argument applies to g=4q3, provided jee and toqg=q, 
1 
| provided p< a If, however, g=q, and pe or g=g¢, and jee then 
2 1 2 
the solution is in the form (51). 


Solutions of Cases (d) and (e) can be obtained from the solutions of Cases (b) 
and (a), respectively, by the transformation mentioned in the preceding section. 
/ When @, is expressed as a function of time, the corresponding expressions for 
@, and @, can be obtained from (8) and (9), respectively. 
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Introduction 
In an optimum control problem we are given a real ordinary differential 
equation system 


a St i PR ee ae ee, P=, De ns yt 
which is a mathematical model of some physical process. The problem of control 
is to select the real functions w(t), 7=1,2,...,m (control variables) on an 


interval of time, fj <¢<#,, such that the solution x (¢) moves in a prescribed 
manner on ¢,<t<Xt#,. The quality of this choice is measured in terms of a 
performance index. For example, it could be required that the «‘(¢) move from 
a prescribed initial point to a prescribed moving target G'(¢) in a minimum 
interval of time by choosing the w(t) from an appropriate class of controls; 
or it could be required that the x'(¢) move to G'(é) in a finite interval of time 
by using w/(¢) in which the energy for control is to be a minimum. Usually the 
performance is measured by a functional which depends on the control variables 
wi (t) and the controlled variables x‘ (t). 

We shall consider the problem of existence of an optimal control. This 
problem has been solved for the case of linear differential equations in [2], [3], 
[9], and for certain nonlinear equations in [8]; however, our treatment includes 
these earlier results. In Section1, Theorems 1 and 2 state conditions such that, 
if there exists one allowable control which does the prescribed task, an optimal 
control will exist. The results of Section 2 are concerned with establishing the 
existence of an allowable control which accomplishes the task for various forms 
of the differential equation system. 


1. Existence of an Optimal Control 
Consider the differential system 


i 
1) = = ft, e... ew ay: MiG WS ota D 

where* f*(¢, x1,..., x", w,..., u”)=f*(t, x, u) together with 
oft 


Seu); th 1,2,.0.50 


* In vector notation *=(x1,..., #”) and |4|= oe 
i=1 
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} are real continuous functions in R! xR” xQ, where R” is the real n-dimensional 
| number space and Q is a non-empty compact subset of R”. 


For each choice of the function 


u(t) =(#(t),...,u"()) on —w<t<t<i<o, 


/ as a measurable vector-valued function with a graph in 2c k”, the differential 
| system 
1) u(t) Raa elt, x, w(t)) i De en 


| has a unique absolutely continuous solution, or response, x(¢) on a subinterval 
of ¢)<¢Sz%,, with a prescribed initial point %»=x(é)). This is the conclusion 
of the CARATHEODORY existence theorem for differential systems [4]. Note that 
u(t) is measurable if and only if every component w!(Z), ..., w(t) is a real-valued 
(Lebesgue) measurable function. The response x(¢) is continuous, and it has a 
derivative, except on a set of measure zero, such that the differential system 
1).«) 18 satisfied almost everywhere. If the finite interval )<¢<z, is degenerate 
so t)=1,, then the response is just the single point x (tj) =x. 


Definition. A control (or steering function) for system 1) with prescribed 
non-empty compact set Qc R” and prescribed initial point x)€R” is a measur- 
able vector-valued function u(#), on a finite interval 4) <¢<4,, with w(t) CQ, 
_ such that the response x(¢) with x(f)= x» is also defined in R” on 4)<t<4,. 

We shall be interested in those controls such that the response x(t) travels 
from the prescribed initial point %»=~%(t)) to a given moving target G(¢). For 
each ¢ on a given finite interval t,<¢<71,, we specify a non-empty compact 
target set G(t)c R”. Moreover, G(t) varies continuously with ¢. Here we use 
the Hausdorff metric distance between two non-empty compact subsets X and 
Y of R” which is the smallest real number d=d(X, Y) such that X lies in the 
d-neighborhood of Y and Y lies in the d-neighborhood of X, cf. [J]. If G(é) 
is a point for each ¢, then the target is a continuous curve. If G(é) is a constant 
compact set, we have the regulator problem where the target is fixed, cf. [10]. 

Let us give /(é, x1,..., x”, w,...,w”) as a-real continuous function on 
R! xR" xQ and define the cost functional of a control w(t) on 4y<¢tSi#,, with 
response x(t), by 


ty 
C(u) = f f(t, x@), u@) at. 
to 
If f(t, x, wu) =1, then C(uv)=t,—%), and the cost of a control is just the time 
duration over which it acts. 


Definition. Given the control problem 


Kia (tx, ys P= Iy Bo ocay ls 


Git) CRA, ON Tess 6S 5 
G 
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as above. Define A=A (f(t, x, 4)... ( %, 4), Q, %, G(é)) as the set of all 
controls w(é), on various subintervals 4;S¢S#, with tShS4S%, such that 
x(t) =% and x(4)€G(4). 
A control #*(¢) in A is called optimal in case 
C(u*) < C(u) 
for every control w(f) in A. 

We now prove the basic existence theorem for optimal controls. The examples 
following the theorem show that an optimal control need not exist if the hypo- 
theses of Theorem 1 are not upheld. 

Theorem 1. Given the control problem 


a) #=fh 2,02 0,.., w= DME a) uh for i=1,...,0 and 
Pw A 


i oh’ 
j=1,...,m with g(t, x), Ajt, x), and oy G+), ge (5:5 RHA pda peer 
tinuous on Rix R", 

b) a non-empty, convex, compact restraint set Oe 

c) the imtial point xyER", 

d) the continuously moving non-empty compact target set G(t) on the finite 
interval T StS, 

é) the cost functional 53 
C(u) = f P(t, x), u(é)) at, 

fo 
where f(t, x, u)= 2° (t, x) +h? (t, x) u’, and g(t, x), and h} (t, x) are continuous on 
Te hae 

Assume the set A of controls with responses traveling from xq to G, as defined 
above, 1s such that: 

A) A is non-empty, 

B) there exists a real bound B<oo for all responses x(t) corresponding to A, 
that is, |x(t)|< B uniformly for all responses. 

Then there exists an optimal control in A. 

Proof. Since A is non-empty and the corresponding responses are uniformly 
bounded, infC(u)=m>-—oo for all w€A. Either A is a finite set, in which 
case the theorem is trivially true, or we can select a sequence of controls 

u(t), on <t=e) from A, 
with C(u™) decreasing monotonically to m. Select a subsequence (without 
changing the notation) such that 
ise and Ma monotonically. 
Let us take the case where 
1 i = Pee hor | Og 


and consider the other cases later. Then each w) (¢) is bounded and measurable 
on the interval fj <¢<¢* and thus belongs to the Hilbert space L,[¢*, t*]. We 


assume tg <¢f; for if 4 =df, then m=0 and x €G (¢*) so every choice of control 
is optimal. 
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A closed ball in Hilbert space is weakly compact [5], and thus we can select 
a further subsequence u(t) with 


| u™ (t) > u*(t) weakly in L,[é, e¥]. 
We next show that 
u*() on t Sissi 
belongs to A. 
Now a compact convex set QR” is precisely the intersection of a finite 
or a countable number of closed half-spaces, cf. [J]. Let 


a,y+b20 


be one such closed half-space in R”. Let M be the subset of [¢g, ff] for which 
| u*(t) lies in a; y'+b<0. 
If M has positive measure, 

tt 


J [a;w*() +0] pn (@) dt<0 


) where Y» 1s the characteristic function of M. But 


ti i 


lim f [a,u®)*(é) + 6] @,, (2) Mees [a; u**(t) + b] —,, (t) dt. 


k-—>oo tt 


_ This is impossible since the left members are each non-negative. Thus M has 
measure zero. Since there are only a countable number of closed half-spaces 
considered, u*(¢) lies in 2 except on a set of measure zero. Redefine w*(¢) on 
the exceptional null set of t¢ <t<#f so that w*(t) lies everywhere in Q. 


Next consider the response for w*(¢) on t{ <¢<tf. The response for 
wu (@) on Masts 


is x) (é). Then (using vector notation) 


lim x) (¢) = «*(0) weakly on [i, 4], 
k—>oo 
lim g (é, x (t)) = g*(¢) weakly on [i,t], 
k—>oo 
lim h(t, x" ()) = h*() weakly on [7, ef], 
k—>oo 
lim h (é, x (t)) w& ) = *(t) weakly on [f, ef]. 


Then, for each fixed ¢ on ({ <t<7¢f, 
ty 


lim x(t) =x + lim f [g(s, x (s)) + h(s, x (s)) we (s)] ds + 


k->oo . R—> 09 4(k) 
v 


+ lim f [g(s, x (s)) + h(s, x (s)) u™(s)] ds. 


k— oo 1* 
0 
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Thus 


k>oo 


lim 2 () = xy + f Lets) + *06)] ds. 


Therefore 
lim x) (¢) = %(0) 


exists, for each fixed ¢. Also %(#) is absolutely continuous on [fg, tf] and % (t}) = %9- 


Using the Lebesgue convergence theorem [6], we find 


so £(t)=«*(t) for almost all ¢ on [¢j, ¢f]. Changing the definition of x*(/) to 
be precisely %(¢), we now show that x*(t) is the response to the control u*(t). 


Now 
x*(0) =a + Jim J [g(s, x (s)) + A(s, 2 (s)) u(s)] ds, 


so 


+ h(s, x*(s)) wu (s) — h(s, x*(s)) u*(s) + h(s, «*(s)) w*(s)] ds. 


Since w'")(s)€Q, which is compact, and since h(s, «\*)(s)) >A(s, x*(s)) almost 
uniformly by EGororr’s theorem [6], we find 


a* (i= Xo ata [g(s, *(s)) + A(s, x*(s)) w*(s)] ds. 


Therefore x*(¢) on ¢{<t<tf¥ is the response to the control w*(t). 
Now 
x 7) EGE) foreach k=41,2,3,.... 
So 
a (Et) = Jim [2 (et) — x (0) + 2 YY], 
and 
x* (tf) = jim, [«®) (e)] ! 
If x*(f*) were not in G(é), then there would exist a neighborhood N of the 
compact set G(ff), so that «*(t*) is not in the closure of N. But G(t) CN for ¢ 
sufficiently near #*. Thus x" (t)€N for large k and yet x*(f) is not in N. 
This is a contradiction, and therefore x*(t*)CG(¢*), and the control w*(f) on 
tg <t<if belongs to A. 
Now compute the cost of w*(t). Here 
10h) 


C(u™) = f [g(é, x (4) + nO (z, x (6) 4 (0)] de 


16h) 
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and 
Ii t* 


lim C(u Nal (t, «*(é)) dt + jim i (¢, x ()) w)4 (0) dz. 


k-—>0o k—> 00 px 


{ Just as above we contains 
jim C(u™) = C(u*) =m. 


i Therefore u*(¢) on ¢¢ <t<¢# is an optimal control. 
We return to the assumption 
i) > 2, = t) +4* monotonically. 
) Suppose we do not have 
but instead, for example, 
ge and ah? t 


and the other cases can be treated similarly. Extend each control wu (¢) (at 
least for all large k) to the interval i¢ <t<¢f by defining uw) (t)=u,, a constant 
vector in Q, for ‘<t<7¥. Again define the weak limit in Q, 
u*(t) = jim it CN On Mtg sb UTA 
— oo 


' We must show that each response 
x(t) on Mts 


can be extended to the interval #’<t<z#* using the extended controls. Then 
we shall show that w*(t) on #§ <t<¢f is in A and is an optimal control. 


| It is easily seen that all the compact sets G(¢) on ty <¢<1, lie within one 
sphere, S(0, 9) of radius g, centered at the origin. Thus, for (¢# — ¢{”) sufficiently 
small, each x(t) lies in S(0, @) and has an extended response on if <t<## 
which lies in S(0, 2@). Also 
jim [a eP) — a (#)| =0. 
Just as above we find that 

lim x (¢) = x*(2) 

k—>oo 
at each ¢ on t* <t<##, and moreover x*(#) is the absolutely continuous response 
to the control u*(z), and x*(t)=%p. 


Now 
(i) €G (AP) 
and 
jim |x* (et) — x*()| =0. 
Thus 


RT) = jim, a*(”) = Jim [2 (e?) — x (AP) + o*(4?)] 
or 
x*(t*) = lim x ¢). 


k—> oo 
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Thus 
x* (tf) EG (4), 


as required. Therefore, w*(¢) on fei ee ne, 
Finally, 
1h) 
C(u™) = f [g9(é, x (t)) + AP (t, & (0) w7(0)] at 


2) 


approaches the limit m as k->oo. As above we compute 


lim C(u™) = C(u*) =m, 


k+>oo 
and thus w*(¢) on ¢¢ <t<¢f is an optimal control. Q.E.D. 


Remarks. Consider A(t)), the subset of A for which the control w(¢) and the 
response x(t) initiate at a fixed ¢,. If A(f)) is non-empty, and if the responses 
x(t) for control in A(¢,)) are uniformly bounded, then there exists a control in 
A (ty) which is optimal relative to A(f)). The same applies to the set A (fp, 4) CA (to) 
where the control time interval f;<¢S<#, is fixed. If the differential system 
and the cost functional integrand are time-independent, each control in A has 
the same response as some control in A(t), after a time translation. 

If C(u) >m>— oo for all controls wu in A, or in A(f)), then one requires only 
the uniform bound |«(f)|<B for responses x(f) corresponding to controls u (i) 
with C(w) near m. 

Also f'(¢, x, w) need only be defined and satisfy the hypotheses of the theorem, 
for StS, «COCR, and uc Q2C R”, where © is an open set in R” which 
contains the initial point %), the moving target G(¢), and all the responses of 4, 
or of A(¢,), in a compact subset. 

The hypothesis A) of the theorem concerns the domain of controllability for 
the problem, as will be discussed in a later section. The hypothesis B) is satis- 
fied if 

[f° &, wh ee; ea ee ee 
or if 


”) 
re (2,26, 20) = ee, the me Foe Sag f 


axk 


for some real «, in [t), %]] xR" x. Thus B) is always satisfied if g(t, x) and 
h;(t, x) are linear in x. 

The following examples illustrate situations where the optimal control fails 
to exist or is not unique. 


Example 1. 
4=sin2uu, Y=cos2au, #=—4 in R’, 


The initial point is (0, 0,1), and the target is the fixed point (0, 0, 0) on the 
time interval O<?<¢,<2. The restraint set Q is —1<u<1. The cost functional 
by 


is C(u)= Jf (x+y?) dt, and we consider the set of controls A(0). 
0 
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| Each control w(¢) in A(0) is defined on 0<¢<1. Consider the controls u (2) 
‘such that 


sin 27 w™ (t) = sin 2a kt, 


cos 2 u™) (t)=cos2akt, for k=1,2,3,.. 


)Such piecewise continuous controls are easily constructed. The corresponding 
responses are 


(k) 1—cos2zkt (k) sin2akht 
4 (t) = pees US 
(4) fee) 


|Thus «“) (1) =0, y (1) =0, 2 (1) =0. The cost for each w(t) is computed to be 


Cu) = f eat 1 


Dane? tise" Deep 


jim Cu) =0, 
and m=O is the infimum for all C(u) with uw in A(0). Yet there is no optimal 
u*(t) on OX¢t<1 for which the cost is 
uf 
C(u*) =f (x*2-+ y*?2) dt=0. 
0 
For such an optimal control w*(t) the response needed is «*(t)=0, y*(é)= 
This implies 
sin 272 u*(t)=0 and cos2zu*(t)=0 


for almost all ¢. But this is impossible, and hence there does not exist an optimal 
‘control for this control problem. We note that the coefficient functions of the 
differential equation are not linear in u, but hypotheses A) and B) of the theorem 
are satisfied. 


Example 2. 


ah a he ee A TL ICS 


The initial point is (0, 0,1), and the target is the fixed point (0, 0, D) on the 
time interval 0O<¢Si, <2. Ties Qas the compact but non-convex circle wj-+u3=1 


in R?, Again take Ci et. (x? y?) dt, and consider the controls 4(0). Using 


0 
the controls u{*)=cos 2akt, uf) =sin 2akt, R=1, 2,..., we find inf C(u)=0 for 
u€A(0). Yet there is no optimal control in A(0) which yields a cost of zero. 


Example 3. 
(Ae, See ee ain ke. 


The initial point is (—1, 0), and the target is the fixed point (1, 0) on the time 
interval O<X¢<t,<2. The restraint set 2 is OSuw<2. The cost functional is 
t 


< 1 
C(u) = f (2—y) dt=f (2—¥) dx. Each control w(t) in A(0) is defined on 0S#S2 
0 —1 


and yields a response x(t) =¢—1, y(é). 
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Every response y(t) satisfies the inequalities 
O<y(x)S—Inx*, for x0. 


Each continuous response joining (—1,0) to (1,0) must lie below the curve 
y= —In x? on some interval. 
1 
Thus C(u)>f (2+Inx2) dx=0 for each w(f)€A(0). But for u(t) =u,=2—€ 
=i 


on 0X<t<2, we compute the response x=¢—1, and 


y (x) =—In — @—9+1]. 


The cost for such a control in A(0) is 
1 


Cw) = f {24m [Fe (2—e)+1]bdx. 


=1 


Thus 

lim Ci) —0. 
Hence inf C(u)=m=0 for w(t)€ A(0). Therefore there does not exist an optimal 
control in A(0). Here we note that hypothesis B) of the theorem does not hold. 


Example 4. 
iy ee ede Pee aN 5 in 


The initial point is (—1, 0), and the target is the fixed point (1, 0) on the time 
interval O<Xt<#,<2. The restraint set Q is —1<u<1. The cost functional 
i 1 


is C(u) = i] as i= { a, . Each control u(¢) in A(0) is defined on 0<¢<2 
0 —1 


and yields a response x (#)=t—1, y/(é). 


Every response y(¢) satisfies the inequalities 


se Ee Ol) = 4S 4= 4. 
The two controls u,(¢)=-+1 and u_(t)=—1 are each optimal and achieve the 
minimal cost. Here an optimal control in A(0) exists, since the hypotheses of 
Theorem 1 are satisfied, but it is not unique. 

We close this section with an existence theorem for Lipschitz continuous 
controls, which is valid even if the control w enters the coefficients f(é, x, u) 
in a nonlinear manner. 


Theorem 2. Given the control problem 

a) #—=fi(t,,..., 0", w,...,u");1=1,2,..., nwhere fi (t, x, u) and a (t, x, w) ; 
1,k=1,2,...,” are continuous in R1xR"xQ, Wee 

b) a non-empty compact restraint set QC R”, 

c) the initial point x)C R", 


d) the continuously moving non-empty compact target set G(t)<R" for each t 
on the finite interval ty<t<7,, 
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e) the cost functional 
by 
=f p(t, «@, 4) at, 
to 


where f(t, x, u) is continuous in Rt xR" xQ. 
For a given positive constant A consider the class A(Lip A) <A of controls u(t), 
each continuous and satisfying a Lipschitz condition 


| u(t) — u(t’) 


jfor all pairs t, t’ on some interval tyStySt, USS. Assume 

A) A(LipA) ts non-empty 

| B) there exists a bound B<oo for all responses x(t) corresponding to controls 
of A(Lip A), that is, |x(t)|<B uniformly for all responses. Then there exists an 
optimal control u*(t)€ A(Lip A), that is, C(u*)<C(u) for all u(t)€ A(Lip A). 


Proof. Assume that A(Lip A) is infinite, and define 


inf C(u) =m>— oo for all u(t)€A(LipA). 


‘Select a sequence u)() on 1) <t¢<7 of control of A(LipA) with C(u™) de- 
‘creasing monotonically towards m as k=1, 2,3,... tends to infinity. Select a 
‘subsequence (still called w' (¢)) with 


i) _»+7* and t-+t* monotonically. 
Consider first the case where 
ay et Rete, Saar, 


and again we omit the trivial subcase where ¢¢ =¢{'. Using Asco.i’s theorem [6], 
select a subsequence of these controls such that 


lim wu (2) = u*(t) 
k—>0o 
uniformly on t¢ <#<f and w*(f) is a continuous function satisfying 
|w*(t) —u*(')| SAlt—?| for oSt, USH. 
We must show that w*(é) on tf StS<éf lies in A(LipA). The graph u*() CQ, 
since 2 is compact. 


Consider the responses x) (¢) of w (é) on ¢) <t<7z, Here in vector notation 
x) (t) — xP) (E ne x*) (s), uw (s)) ds — ic f(s, x'?)(s), w)(s)) ds 
(i) 4) 


for (* <t<t# with k and p>k positive integers. Then 


x") (t) — x) (@)| < Sits, x) (s), wl (s))| ds + f | f(s, x? (s), w®(s))| ds + 
) 
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Take k so large that |#% —7| and | ¢* —¢f| are very small and | u (?) — '?) (6)| 
is uniformly small on #<t<#*. Use the uniform continuity of f(¢, x, ~) in 
TSt<t, |x| SB, wEQ, and compute 


t 
|x () — x) (| S = a [Its x) (s), a (s)) — f(s, x™(s), w®)(s))| ds + 
tt 
+f LAs, (5), (5)) — Fs, x9), lP(9)] ds 


ty 
so 


Therefore CAucHyY’s criterion yields 


lim x) (¢) = x**() uniformlyon #*<is #éF, 


k—>oo 
and x*(¢) is continuous on ff St<¢F. 
Now 


IO Ve eae li (s, x (s), u™ (s)) ds + se: (s, x (s), 0 (s)) ds 


10h) 1 
on if =t= 4 Thus 
t 
x*(t) = x9 + f f(s, x*(s), w*(s)) ds, 
iy 


and x*(/) is the response to the control w*(t) on if <¢<##¥. Clearly x*(t*) =x 
Also 


0- 


x(t) = lim x (4) = lim [xe @) — x @) + 2 (e8)] 


k—> oo k—0o 
SO 
a* (ef) lim a). 
—>00) 
Since 
x (AP) EG (), 
we have 


a* (it) €G (tf). 
Thus w*(¢) on t¢ <t<¢* belongs to A(Lip.A). Furthermore 
C(u*) = hm Cu) =m, 


and u*(t) on tj <t<## is optimal in A(Lip A). 
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We return to the assumption that #{") >¢# and ¢\—2*, monotonically. Suppose 
we have tf? <#¥ and t\) <¢*, for example. Extend each control w™ (t) to) <t<t 
py defining wu () =u (¢) on tt? <t<#*. For large k we can also extend the 
s:orresponding responses x") (¢) to the interval #)<#<¢*. Again 

lim u™ (¢) = u*(t) uniformly on i* <t<#F, 


k—+>oo 


jand u*(t) satisfies a Lipschitz condition with constant A. As above we show 
| jim x) (t) = x*(t) uniformly on #* <t< ¢¥ 

jand that x*(¢) is the response to u*(t). Also x*(tj) =x, and x*(t*) CG (t¥). Thus 
ju*(t) on tg <t<t## belongs to A(Lip A). Again we compute that w*(¢) on ¢*{ <t<t* 
jis optimal in A(LipA). Q.E.D. 


Remark. The remarks following Theorem 1 are also applicable here. 


Corollary. Assume the hypotheses of Theorem 2 for each Lipschitz constant 
p=1, 2,3,..., and let u(t) on tf) <t<t? be optimal in A(Lip p). In addition 
assume 

1) the target G(t) on ty StS, ts a fixed point G, 

2) for each s>0 there exists a measurable control u,(t) on t<t<t) in A such 
that C(u,)Sm-+s where m=intC(u) for uEA and also t!) <1, —« for some e>0, 
independent of s, 

3) there is a uniform bound |x (t)|< B< oo for all responses x(t) corresponding 
to measurable controls in A, 

4) for each s>0 there exists a neighborhood N of G such that each point in N 
can be steered to G by a C! control having a prescribed imtial instant th<t,—€, 
a duration <e, and a cost <s (this condition is considered in Theorem 4). 


“a lim C(ut) =m=infC(u) for wed. 
p— co 


Proof. It is clear that 
| Gi ie Cena P= 1,2,3; 0c 


. 
so lim C(u}) exists and is not less than m. But given a measurable control w, (¢) 
poco 


in A, there exists a C® control which approximates wu, (t) almost uniformly so 
that the corresponding response has an end-point in N. But then x» can be 
steered to the target point G by a control w,(t) in A(LipA) for some A>0, 
and with C(w,) arbitrarily near to m. For p>A we thus have 


C(u4) 2 C(up), 


so lim C(us)=m, as required. Q.E.D. 


p— co 
2. Domain of Controllability 
In this section we investigate the nature of hypotheses A) of Theorem 1. For 
simplicity we consider the problem of steering an initial point to a fixed target 


point, say the origin. We still maintain that the restraint set Q is compact, 
but in this section we require that the zero vector is in the interior of Q. 


48 E. B. LEE & L. Markus: 


Definition. The domain of controllability @ for the differential system 


for a given restraint set Q¢ R”, consists of the set of all points x),€ R” for which 
there exists a measurable control u(t) (2, defined on some finite interval, steer- 


ing x) to the origin® of R”. 


Theorem 3. Consider 


ea Ph, cia We | aa Aan te 


where** fi(x, u)EC! in R"xQ. Assume the vector u=0 1s an interior pornt of the 
compact restraint set QCR”. Assume 


(0,0) =O » 4240... 


If the matrix a (0,0) has rank n, then the domain of controllability @ is an 
U 


open connected subset of R” containing the origin. 

Proof. If the origin of R” is an interior point of @, then it follows from general 
continuity arguments that @ is open and connected. By continuity there exists 
a neighborhood N of the origin in 2 such that for any measurable vector 1 (¢) 
on —1<¢<1 in NN the corresponding solution x(¢) initiating at %»=0 will be 
defined in R” for —1<$#<1. 

Let 

GE wy aka): fone Meta i=1,...,n, 


be the solution of the given differential system for constant (w#,..., w”) in N 


and with gy (0, v)=0. We shall prove that — (ty, 0), for some #)< 0, has rank n. 
Hence the differentiable map ao 
U> x= (tb, 4) 


carries a neighborhood of #=0 onto a neighborhood of x»=0. Then there is 
a neighborhood of x)=0 which can be steered to the origin in time —f,>0, 
and each required control is just a constant. 


Now 


oY =F (p(t,u),u), g(t, 0) =0. 
We compute 


a (2%) = of 


ot\aw] — axk pes 


For 4%)=0 write the matrix 


z() =" (0). 


* This is sometimes referred to as null controllability. 


** That is, f(%,u)€C® in R”xQ and has a C! extension to R"x O, where Q is 
an open set containing Q. 
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dz 


| 
| 
ap = 1(0,0)2+4,(0,0), (0) =0. 


| Then 
. 


| Then 


For a small |é,| 


i has rank n, since /,,(0, 0) has rank 2. Q.E.D. 


Corollary. Consider a linear differential system with real constant coefficients 
ti = At a* + Bil Ula, ede Ha, oo, ee 


| Assume the compact restraint set Q 0 R” contains the zero vector in its interior. 
\ If the matrix B has rank n, the domain of controllability €@ is an open connected 
» set in R” containing the origin. 
A more difficult problem arises when the control has fewer degrees of freedom 
, than the response, 7.e. when m<m. The next theorem considers this case. 
Theorem 4. Consider 

Ret ween OU cca, a0") (eer 


where f(x, u)ECCt in R”xQ. Assume the vector u=0 is an interior point of the 
compact restraint set 2Q°R”. Assume 


1) f'(0, 0) =0 (Ne he arse) 


and 
2) there exists a vector vER” such that Bv les in no invariant subspace of A 
with dimension <n —1, where peso (0,0) and A = (0,0) are real matrices. 


Then the domain of controllability € 1s an open connected subset of R" contain- 
ing the origin. 

Proof. Consider a neighborhood N <Q of the origin in R” and t,>0 such that 
for each measurable vector w(t) on —1,S¢<71, with w(t) CN the corresponding 
response x(t) with x(0)=0 is defined on —1,S/S 7}. . 

We shall define a family of controls in N 


ul (6 E) == UE, Eyes Yen) ASS Gt 


on |¢|<7, and |é|<e for some e>0. Take 


u(t) =0[XGH)h4+X(6S)a+--+X(,2)8,| 


where v is the vector mentioned in hypothesis 2), 


4 on |t| Sh 
OCT = 
if) {5 on |t| >A, 
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and ¢, is a small number with 0<|t,|<1%. For convenience we choose this piece- 
wise continuous function X(t, h), but it would be easy to use a C® function 
which almost uniformly approximates X(é, h), for each fixed h. 


Note w(t, 0) =0 and also 


2 =0X(t, 3), f= 4 le worsen 1. 
ag 


is differentiable except at ¢=1,/k. 


Let 
MD (b-Gaaee val, ls (= lear 


be the response corresponding to w(#, €) with 
a (0; eye Os tA, cna, Me 


Note. that# (4,0) =0, t=4,.-.; %, simee f'(0, 0) =0. Now x(é, €) is continuous 
in (n+1) variables. Also, for each fixed ¢, the map 


& > x(t, &) 
is a differentiable map of a neighborhood of the origin of R” into Rk” with the 


origin fixed. If 


: 4 
dlis a (£0) a tke Atos 0 
Ss 


is non-singular, for some fixed t<0, then the domain of controllability @ is 
an open connected subset of R” containing the origin, as required. 


Now, in vector notation, 


x(t, €) =f Hx £), u(s, €)) ds, 


and so 


ra) } 
g¢%O=[|AZ6O+B he (s, 0)| ds 


or 
t 


z(t) =f [Az(s) + Bug(s, 0)] ds, 


0 
so z(t) is continuous near {=0. Also dz/dt exists, except at t=1,/k, R=1,2,...,%. 
The k™ column of the matrix z(t) satisfies 
t 
tay) = [ |Azay(s) + Box, rata 
0 


or A=41; 2, 2:5, where the subscript k=1, 2,...,” designates the column. 
Thus we have 
da(n 
dt 


AG. t Bo X(t, 2) | es Ieee 
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t 
a vA SA! ue 
Zn (t) = € fe BoX(x, hj ds 


0 
or ty/k 


raea= tiie + Bods: 
0 


a %,)(0)=0. But the solution of this linear differential system is 


Thus, for |4|<|t]<y, 


t 3 4 

way) =e4)hI— A+ 2 At 2 aay | Be, 
t t 2 3 

wal) e427 — GE 4 4 (AIA 42 ...] By, 


hy (t,/”)? A4 (t,/n)8 A2 | Bo 
: 31 | 


We shall show the linear independence of the vectors (e~'4) zq) (4), ..., (e€7'4) 2, 
The Vandermonde determinant 
ty MOA 7 


t[2 (6/2)... (6/2)"| 4. 6, 


tn (ty/n)® ... (t/n)" 


Thus it is sufficient to show the linear independence of the vectors 


Bv+0O(t) 
ABv+0(4) 
A? Bv +0 (t) 

A" By + O(t), 


where lim O (4) =0. Since the vectors Bv, A Bv,..., A” *Bv are linearly in- 
t— 


| dependent, we conclude that the vectors 2) (4), 22) (t), ---, 2%) (¢) are independent 
and that the matrix z(#) is non-singular. Q.E.D. 


Corollary. Consider 
Ree Yala eels a ah hey teen. 1 

where f(x, u)€Ct in R” xQ and the origin of R” ts an interior point of the compact 
set Q. Assume 

te 0-0) = One tes ie I 

2) there exists a vector vCR” such that Bu les in no invariant subspace of A 
with dimension <n—1. 

Assume also that 

Reif (Ce wer Ore 92,0) t=1,...,” 

is globally asymptotically stable* towards the origin of R". Then the domain of 
controllability @ is the entire space R". 


* See reference [12]. 
4* 
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Corollary. Consider 
= , —1 _— 
GOL ge (HX, 5 ce) YE hae, ao oad Be cag 


where the coefficients a;(x,, X2,-.-, %,)€C1 im R" and |u|Se, for some e>0. 
We write the corresponding first order system 


Xy =X, 
S) bee 
<= Ay —4(%, COL Kn) Lo ret? aaa Ay (%4, Plas Xn) My U. 


If the system 


Ay — Xa) 


Hace A 555 Ne drab tg) Mat a ae Tig Fates x By) xy 


is globally asymptotically stable towards the origin, then for S) hypotheses 1) and 2) 
of the theorem automatically hold and the domain of controllability for S) as all R". 


Corollary. Consider 
TH Sl PRISE ES oon oH 4:42 oh 


where f(x, u)€C! for xER" and uCUCR™. Assume that for each point x CR” 
there exists an interior point u(x%») of U such that 


DACA ONC = Or VA oueeeg Ht 
and 


2) there exists a vector ve R™ such that Bu lies in no invariant subspace of 


A with dimension <n—1, where B= - (% 5, (%»)) and A = s. (%, U(X) are 


uw x 
veal matrices. Then for each pair of points x, and xy of R” there exists a piecewise 
continuous controller u(t) on t<t<t, in U which steers the response x(t) from 


(Gi) — Ay 00H (tn) = eo 


Proof. In the proof of Theorem 4 we note that there is a neighborhood of the 
origin which consists of points that can be steered to the origin by a piecewise 
continuous controller in a finite time. Also there is a neighborhood of the origin 
consisting of points to which the origin can be steered by a piecewise continuous 
controller in a finite time. 


The same properties hold for each point % in R”, after translating both x 
and uw as required. 


Now let S be the subset of R” consisting of points to which x, can be steered 
by a piecewise continuous controller in a finite time interval. Clearly S is both 
open and closed in k”. Thus S=R” and x, can be steered to x5, as required, 
Chery (manly 


_ Remarks. In the previous corollaries we can replace the hypothesis that 
X'=f'(x,..., x", u,..., vw") be asymptotically stable for w/=0, 7=1,..., m, by 
the hypothesis that #'=/*(x!,...,%”, #,..., wu”) be asymptotically stable for 
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|some appropriate choice of u/=w (x1, ..., x"). Asa particular choice of «!=w (x) 
to satisfy the condition of asymptotic stability consider a Lyapunov function 
'V(x)20, V(x)E€C, with the choice of #€Q2 which minimizes Es f' (x, u) 
lf this minimum is negative definite in 0<|x|<oco for some V(x)>0, with 
V(x)—>oo as |x%|—>co, then the domain of controllability @ is all of R”. 


We next consider the domain of controllability for certain second order 
: differential equations 


x +(x, 4%, u) =0 


which we write as the first order system 


EAR Nah EN 
in the phase plane. 


Theorem 5. Consider the differential equation 
% +f (x,.%, 4) =0 
» with f(x, y, uw) €C! for all (x, y) and wE€Q, where Q is a compact interval containing 
zero as an interior point. Assume 
ea) (Os0;,0)s==10 
b) f(x, ¥, 0) >0 and f,(x, y, 0) >0 everywhere 
Cy fMONOKO) 20. 
Then the domain of controllability @ 1s the entire (x, y)-plane. 
Proof. The system 
fay 
y =—f(*, ¥, 0) 
O 1 
a, ps eo if : 
have negative real parts for all (x,y), the system is globally asymptotically 


stable to the origin, cf. [12]. However the solutions will not reach the origin 
in a finite time. 


has the variable Jacobian matrix J =|( ) Since the eigenvalues of / 
4 


Now 


y =|( 0 1 
—f,(0, 0, 0) Saif AO; 0, 0) 
and 
0) 


| caer a 0, ») 


Take v=1 as a 1-vector; then 


gels 0 ==) 10; 0, 0) 
Ze a 0, 0) ta ty 


are linearly independent. Thus the domain of controllability to the origin is 
an open set. Therefore @ is the entire (x, y)-plane. Q.E.D. 


and ABy=( 
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Corollary. Consider % + f(%) + g(x) =u 


where f(y) and g(x) €C! for all (x, y) and |u|Se, for some e>0. Assume 


a) ¢(0) =g(0) =0 

b) f(y) >0 and g'(x)>0 everywhere. 
Then the domain of controllability @ is the entire (x, y)-plane. 

We now turn from the asymptotically stable case, where the domain of 
controllability @ is the entire space, to the conditionally unstable saddle-point 
case, where @ is a strip in the plane. 


Lemma 1. Consider 
QZ) % + f(%) — g(x) =0 


with f(y), g(x)EC! and f(0)=g(0)=0, f(y) >0, g'(x)>0 for all (x, y) in Re 
Then the solution curve family of G) in the (x, y=) phase plane 1s topologically 
equivalent to the solution curve family of the linear equation 


%+%—*K=0. 


Proof. Write the equation ZY) in the phase plane as 


~) ey 

PSE Sy) 
We seek a homeomorphism of R? onto R? which carries the solution curves 
(sensed but not parametrized) of ) onto those of 


L=y, 
L) teat 
OF eee 
In Y) the origin is the unique critical point. The variational equation o 
SF) at the origin is 


veil 
y=s'(O) xf) y, 


with real eigenvalues of opposite signs. Thus ) and ¥) are topologically equi- 
valent near the origin. 

Thus ) has exactly four solution curves which have the origin as a limit 
point and these approach the origin with the same directions as for ¥), that is, 
one of these solutions of ) lies, near the origin, in each of the four quadrants. 
Call these solutions of ), which will be shown to be separatrices in the sense 
of Markus [13], by the numerals I, II, III, IV corresponding to the quadrants 
in which they approach the origin. By examining .Y) on each of the coordinate 
axes it is easy to see that each of I, II, III, IV lies entirely in the corresponding 
quadrant. 


Along I we have dy g(x) —f(y) 
Ox v4 


which is bounded from above on each compact x-interval. Hence I is single- 
valued over the entire positive x-axis. A similar argument shows that III is 
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_single- valued over the entire negative x-axis. The separatrices II and IV are 
each single-valued over a segment of the x-axis, and they extend to infinity, 
j that is, x24 y2—>0o on II and IV. 

_ Consider a solution curve of SF) in the plane sector bounded by separatrices 
/ILand I. Here ~=y>0 and, in the second quadrant y<0. Thus each solution 
curve in this sector must intersect the first quadrant and must exist in the first 
quadrant for all larger values of x to the right of any such intersection point. 
) Hence the solution curves in the sector between II and I are linearly ordered 
| (using the inclusion relation for the regions above the curves), and they form 
a parallel family, as described by Markus [13]. A similar analysis holds for the 
sector between III and IV. 

In the sector bounded by IV and I each solution of Y) intersects the positive 
| x-axis, which thereby serves as a transversal for the solutions in this sector. 
A similar analysis holds in the sector bounded by II and III. 
| Thus the only separatrices of .) are the critical point at the origin and 
‘ the solution curves I, II, III, and IV. Therefore by the general theory of separatrix 
configurations [13] we find that ) is homeomorphic with Y), as required. Q.E.D. 

Remark. Call the solution curves II and IV, leading towards the critical 
point of ) with negative slope, the principal separatrices, and call I and III 
the minor separatrices of .). The principal separatrices of .Y), together with 
'the critical point at the origin, form a topological image of a line in R? which 
‘is single-valued over the entire y-axis. The same holds for the minor separatrices 
of /) over the entire x-axis. 


Lemma 2. Consider the two differential systems 


aes 
Ve) I 19) 


(x)=g(x)+C,, for positive constants C,, are in C1 in R*. 
Assume {(0)=g(0)=0, f’(y)>0, g'(x)>0 everywhere in R?, and |g(x)|—>0o as 
aloo. 
Then each of ¥,) and F.) is topologically equivalent to the linear system, 
in the phase plane (x, y=%) 


gL) eae 
= RVG 


where f(y) and g, 


The principal separatrices Il, and IV,, together with the critical pornt, of 
S,) form the topological image of a line which separates the plane im two. At each 
fixed ordinate the principal separatrix or critical point of F_) lves to the right of 
the corresponding principal separatrix or critical point of F,). The open band 
B between the principal separatrices and critical points of F,) and SF.) 1s homeo- 
morphic to an open band between parallel straight lines in the plane. Also each 
minor separatrix of S,) intersects the principal separatrices of F,) in at most 
one point, respectively (see Figure 1). 


Proof. By Lemma 1 the topological configuration of the solutions of %,) or 
of £) is that of Y). It is clear that the critical point of %) lies to the right of 
the critical point of Y,), on the x-axis. By the preceeding remark the principal 
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separatrices II_ and IV_, together with the critical point of “), form the topo- 
logical image of a line which is a single-valued covering of the y-axis. The same 
holds for Y,). 

Suppose II, and II_ intersect. At the intersection point which is furthest 
to the right the slope of II_ is less than that of II,, and this contradicts the 
supposition of an intersection. Similarly IV, and IV_ cannot intersect. 

Therefore the open band B bounded by the principal separatrices and critical 
points of Y,) and “%) is homeomorphic to the open band between two parallel 
straight lines in the plane. In fact, the homeomorphism can be extended to the 
closed bands. 


Fig. 1 


Now on a minor separatrix, say L,, we have ¥>0 and y>0. Hence I, intersects 
II_ in exactly one point. Similarly III_ intersects IV, in exactly one point. Q.E.D. 


Theorem 6. Consider 
Z) E+ H(%) — g(a) =u, 


where f(y), g(x)€C for all (x, y) in R? and —C_SusC, ts the finite real interval 
Q with —C_<0<C,. Assume }{(0)=g(0)=0, f’(y)>0, g’(x)>0 im R? and 
|g(x)|>00 as |x|—>0o. Then the domain of controllability @ in the phase plane 
(x, y=) ts precisely the open topological band B bounded by the principal separa- 
trices, and critical points, of the two systems 


Sf) ee 
y= a(2) £(C.) 0). 


Proof. Consider the open quadrilateral Q in R® bounded by I,, IV, and 
hy The origin 0 lies in Q, and the solution ,) through 0 must pass 
through a point 0, in Q very near III_ before it comes to 0. 


Now take an initial point # in Q. Steer # by the solution of “) until it 
almost reaches the intersection of this solution with II_. Then follow the solution 
of “) around the inside of the boundary of Q to the point 0,. Then switch to 
the solution of %,) which steers the point to 0. 


It p lies in the band B but above the quadrilateral Q, start out along a 
solution of -,) until one almost reaches’ the intersection with II_. Then switch 


to a solution of %) and follow this until the point 0, is again reached. Then 
follow the solution of Y,) into 0. 
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If # es in the band B but below the quadrilateral Q, start out along a solution 
‘of £) until one almost reaches IV,. Then follow “) into Q, and proceed as 
above. 


Now let Z be a point of R*—B. If Z lies in y>0 and to the right of the 
critical point of Y), then Z cannot be steered out of this quadrant by a controller 
)—C_Su/(t)<C, and thus Z cannot be steered to 0. Now Z follows a solution 
‘curve of 


eae 
j=e.(x) fy) +20) 


twhere «¢(t)=wu(t)+C_=20 is measurable on some finite time interval. Thus if 
1Z lies in the open sector bounded by II_ and L_, then Z must lie above a solution 
icurve of “), and hence, using the parallel structure of the solution curve family 
in this sector, we see that Z cannot be steered to 0. If Z lies in the open sector 
‘bounded by IV_ and I_, then Z must enter the half-plane y>0 to the right 
‘of the critical point of %), whatever the controller w(t). Therefore, if Z lies to 
‘the right of the closed band B, then Z cannot be steered to the origin 0 by a 
smeasurable controller u(é) with —C_Su(t)<C,. A similar argument shows 
that if Z lies to the right of the closed band B, it cannot be steered to the origin 0. 

By Theorem 4 the domain of controllability @ is an open plane set. There- 
More 6 — B. Q.E.D. 

Remark. Each initial point in @=B can be steered to the origin 0 using only 
the solutions of Y%.) and “). It is necessary to use only two switches between 
these two systems for each point in @ (even though the discussion in the theorem 
uses three switches in certain cases). 


As a final example consider the van der Pol equation 


for a positive constant yw. In the phase plane this is the system 


AN) 
y=e1— x) y— 2. 


The origin 0 is the unique critical point and is an unstable focus, or node. There 
is a unique periodic solutions S, which is orbitally asymptotically stable and 
which lies between the abscissas x= +d, (for estimates of d,, see [14]). 


Theorem 7. Consider the differential system 


i, 
q ere lynsey a 


for a fixed 4>0 and measurable controls u(t) defined on finite intervals with 
|u(t)|<e, for some e>0. Then the domain of controllability @ to the origin 0 1s 
all R* whenever 

e>d,/2. 
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Proof. For each choice of a constant w on |u|<e the corresponding system S) 
has a unique critical point 0,, which is unstable, and a unique periodic solution, 
which is orbitally asymptotically stable and which lies between the abscissas 


x=u-td,. 
If fie constant u satisfies |w|<d,, then it is easy to see that 0, belongs 
to @. If |u|<d, and if @ intersects the periodic solution S,, then @= k®. 
Now if e>d,/2, we can choose a constant w so |u|<e and yet |2u|>d,. Then 
@ intersects S,, and hence @=R?. Q.E.D. 


be? 
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Finite Dimensionality 
ot the Schrodinger Operator Bottom 


F. H. BROWNELL 


Communicated by M. M. SCHIFFER 


1. Introduction 


We consider in the complex Hilbert space L,(E,,), E,, being n-dimensional 
‘Euclidean space with n=1, the Schrédinger partial differential operator H. 


0 


(1.1) [Hy u] (x) = — V2? u(x) + V(x) u(x) 


over CE, for w©Y, the domain Y of H, being the set of all complex-valued 
‘functions uw over E,, such that VwEL,(E,,), that all partial derivatives of u of 
order <2 exist and are continuous throughout F,,, and that uw and all such 
partials remain bounded over E,, after multiplication by (1+ |a|”)+ exp ($ |a|?) 
4 ae 
“ee ee oa is the Lapla- 
cian, and V is a fixed, measurable, real-valued function over E,, subject to 
‘conditions stated in detail in the next section. These conditions suffice to specify 
a self-adjoint extension H of H,. Let E be the spectral measure of H, and 
let N(A) be the dimension of the (—oe, A] eigenspace E ((— ov, A]) L,(E,,) for 
real 4. Our main result is the proof (under our conditions on V, which imply 
that im (ess inf V(a))=c) that N(A)SNj<+oo over all A<c. By adding a 


+oo |xl|2r 

‘suitable peacuat we may and hereafter do consider only the c=O case for 
simplicity. In the next section we state precisely the conditions on V and our 
resulting main theorem. We remark that these conditions are implied by the 
more stringent ones of IkEBE [6], for which our conclusion here thus holds. 
Actually the variational machinery developed in a previous paper on asymptotics 
[3] enables us readily to reduce our theorem to some well known lemmas about 
ordinary differential equations (also stated in the next section). The stated 
result for (1.1) from these known tools does seem to be new. 


for some m=0 depending on w. Here |x| = De |x; |? 
7 


2. Hypotheses and Main Statements 
Let y,, denote n-dimensional Lebesgue measure for E,,, and let “‘measurable’’ 
mean measurable with respect to the o-algebra of E,, subsets which is the domain 
of 4. With this notation our hypotheses about V are stated as follows. L,(D) 


This research was supported financially by the National Science Foundation, 
grant NSF-G 11097. 
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for measurable DCE, means the complex L,-space formed with respect to ,- 
integration over D. 


Condition 2.A. Let V be a real-valued, measurable function over E,, which 
is essentially bounded over {a||a|= Ry} for some real Ry>O and which for some 
veal y=1 having y>tn satisfies 


(2.4) Sf |V(x)|’ du, (#2) < + 0. 
jx|SR, 

Condition 2.B. Let V satisfy Condition 2.A, and in addition let V (a) = — g (|x|) 
almost (1,,) everywhere over \a|=Ry, where g(r) 1s defined and real-valued and 
possesses a continuous derivative over real r2R, with g(r)=O0 there such that, 
for some finite M and RZ R,: 


(2.2) g(r) SM rr over r€[Ry, +09); 


(2.3) rfg(s)dsSt over r€[R,+oco) if |n—2|=1; 
(2.3)’ hetdstel tial ae over [R, +c) 4] 42 


Notice that both (2.2) and (2.3) are implied if 1>lmsup4r? g(r); likewise 
both (2.2) and (2.3)’ are implied if 1> im sup 4r? ne r)? g(r). Also (2.3) is imphed 
he 


by frg(r)dr<-+ow, and (2.3) is implied by fr (Inr)? g(r) dr<-+tos, since the 
Ro 
right side of (2.3)’ has the asymptotic form 4 fines )?{14+0((nr)+)} as r++ oo. 


Assuming only Condition 2.A for this purpose, we now specify our self- 
adjoint extension H of H, in L,(E,,) as follows. Define H=G,'—w?IJ as in 
Theorem T.5 of [2] (see T.3 and T.4 there also, pp. 959—961), H being self- 
adjoint since G,, is a bounded Hermitian integral operator on L,(E,,) with Green 
function kernel. Also HjCH as desired for extension ((2], Theorem T.6, p. 964 


except for the trivially justified replacement of 2, there by D of (1.1) here). 
Other than these references to [2], we here give no further details of this con- 
struction of H, since for y=2 in Condition 2.A, as is required if dimension n= 4, 
this construction will now be seen unnecessary to specify H. 

Now assume Condition 2.A with y=2 also. Then the condition VwE€L, in 
the definition of D following (1.1) may be dropped as redundant, and &, is 
dense in L,(£,,) and Hy is symmetric. Also under Condition 2.A with y=2 
there exists a unique self-adjoint extension H’ of Hy (by [2], Theorem T.1, p. 954 
with H,CH CH’ and H’ the unique self-adjoint extension of H, and thus also 
of H,). Thus the self-adjoint extension H constructed above must coincide with 
H', and H here may be specified as this unique self-adjoint extension of Veh, 

Thus we have specified our self-adjoint extension H of Hy, using Condition 2.A 
only. With N(A) defined in terms of the spectral measure of H as in the Intro- 
duction, we may state our main theorem as follows. 


Theorem 2.C. Let V satisfy both Conditions 2.A and 2.B. Then N(A)<N, 
over real A<O for some positive integer Ny< oo. 
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By using the variational machinery in [3], we can reduce the proof of 
Theorem 2.C to some well known lemmas about ordinary differential equations, 
chis reduction process making no use of either (2.2) or (2.3) except for the con- 
dition that 0= lim g(r), as is obviously implied by (2.2). In Section 3 we shall 


r—>+00 


describe this reduction procedure and in Section 4 discusss and give the references 


‘or these lemmas from which Conditions 2.3 and 2.3’ arose. 
| For these lemmas we consider for real R the differential equation over (R, + oo) 


a — ul'(t) = f(r) u(r) =— 08 u(r), 


eer by a solution u of (2.4) over (R, + cc) we mean a complex-valued function 
there such that w’’(r) exists satisfying (2.4) at each real r>R. We shall be 
Biferested also in solutions w of (2.4) which satisfy at R a boundary condition 
of the form 


(2.5) 0 = (cos #) u’(R*) — (sin 8) u(R*) 


for a given real #, where the right-hand limits w’(R*) and w(R*) are assumed 
to exist and be finite. 


; Lemma 2.D. Let the function f be continuous and real-valued over |R, +0) 
bad satisfy f(r) =0 there for some real R,>0. In addition let f satisfy either 


(2.6) rff(s)ds<4 


over r= R, or else satisfy f(r) =1r-24 f(r) with 

i(2.6)’ r{f(s)ds<4rfs(Ins)2ds 

lover r=R, (the Lebesgue integral on the left of (2.6)' existing since f(s) ea 
‘Then there exists some real R'=R, such that for every real 0 and veal R= R’ there 
exists at most a single real w>0 having a solution u of (2.4) and (2.5) for R which 


also satisfies 0< f |u(r)|2dr<+oo. 
R 


Finally we record the following well known and trivial lemma (differing only 
in boundary condition from Theorem 3.9, [3]). Solutions are defined as for (2.4). 


Lemma 2.E. Let q be a real-valued, bounded, continuous function over (R, + 0%) 
for some real R, and let the real number A satisfy 


(2.7) 4 < inf g(r). 


r>R 


ien for each real k=O there exists no real-valued solution u of 
(2.8) —u''(r) + q¢(r) u(r) =Au(r) 


over (R, +00) which also satisfies 0< J |u|®diy<-+ co and has finite right-hand 
limits u'(R*) and u(R*) satisfying * 


(2.9) 0 =w'(R*) — ku(R*). 
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In closing this section, we remark that Theorem 2.C taken together with 
our asymptotic formula in [3] indicates, by comparison of (2.2) here with (2.1) 
o (2.6) in [3], that O(r-*) at oo for the potential function is roughly speaking 
the borderline between N(A) bounded or unbounded as A—-0~. In addition to 
this description of the (— ov, 0) part of the spectrum of H, theorems describing 
the [0, --co) part of the spectrum are also known ((Z], Theorem 6.1, p. 573, 
Theorem 6.7, p. 586, Theorem 6.8, p. 591; [6], Theorem 6, p. 20; [7], Theorem 14 
and 1a, pp. 405—406). 

3. Reduction Proof 

In this section we shall prove Theorem 2.C, assuming valid Lemmas 2.D 
and 2.E. These lemmas will be proved in the next section. Thus here we assume 
Conditions 2.A and 2.B to be satisfied; actually, as remarked in Section 2, until 
treating the ordinary differential equations at the end of this proof, we shall 
make no use of (2.2) or (2.3) here except for the fact that _him g(r) =0, obviously 
implied by (2.2). Re 

Now to employ the variational machinery ([3], Section 3 plus Addendum) 
needed here, we consider Dj ={xCE,|b<|a|} and Dj; ={ax€E,|b>|x|} for 
bE(0, +cc), and with D=D or D=E,, we consider QpCI,(D) defined as 
follows. Qp is the set of all w€L,(D) such that u possesses throughout D con- 
tinuous partials of all orders <2, that all such partials possess continuous 
extensions over the closure D of D, and (if D is unbounded) that wu and all such 
partials are each in absolute value <M(1 + |a|*) exp (— 27|a|?) over «€D for 
some M and « in (0, +00) depending on w. Finally, with the function q(x) 
defined either by ¢(x)=V(a) or g(#)=—g({a|) in Conditions 2.A and 2.B, we 
define for such D and q and for each integer p= 1 


G1) t(D.) =sup elt {|W (e)|* + ga) | (e)|2} de, (w)], 


2 
where |V u(x Ss , where the outside supremum is taken over all 


i= 
Sets 5 CLAD) nea es S=p—1 (and thus S=Q, the empty set, if =1), 
and where the inside infimum is taken over all w€Qp such that fw@du,=0 
D 


for each wES and f|u|?du,=1. Finally, define the non-negative integer 
D 


m(A; D,q), for real A and D and q as in (3.1), as the number of integers p>1 
such that t,(D, g)SA. 

With this notation, Conditions 2.A and 2.B together imply for b—1=R, 
the following inequalities over 4<0: 


(3.2) N(A) =m(A; E,,, V) < + 00, 

(3.3) m(A; E,,V) <m(A; Dy, V) + m(A; Dj, V), 
(3.4) m(A; Dy, V) S m(0; Dz, V), 

(3.5) m(0; Dy, V)< +, 

(3.6) m(A; Dy, V) Sm(A; Db, qe), 


where q(x) = —g(|a|) is defined over #€D}. 
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Here (3.2) is proved in [2], Theorems T.7 and T.13; also (3.3) and (3.5) are 
sroved in [3] (see Addendum, second paragraph), Lemma 3.2 and Lemma 3.7 
here respectively, the last using boundedness of V(a) over b>|a|>b—1=R,. 
3.4) is obvious from the definition of m(A; D, q), and V(x) = —g(|a|) =q (a) 
umost (w,,) everywhere over |a|>b=R,+1 by Condition 2.B yields t,(Dj, V)= 
‘p (D3, J) in (3.1) and hence yields (3.6). 

_ In view of (3.2)—(3.6), we see that the desired conclusion of Theorem 2.C 
lollows if we can show m/(A; Dj, qo) bounded over A<0. Since Condition 2.B 


assures that q, (a) =—g (|ar|) possesses continuous first partials throughout |a| > Ry 
md that lim g(r)=0, it follows ((3], Lemmas 3.5 and 3.6 and Addendum, 
second paragraph) for b—1=Ry, g(x) =—g(|a#|) that over A<O 

3.7) m(A; Dg, qo) SNS (Ai Go) < +, 


actually m=N here). Here Nj (A; qo) is defined (see [3] preceding Lemma 3.3 
here) as the dimension of the complex linear manifold (necessarily closed since 
\‘inite-dimensional) in L,(D;) generated by all real-valued wE€L,(Dj;) having 


continuous second partials throughout D;, satisfying 


3.8) — V2 u(a) — g(\al) u(x) = A u(x) 


throughout a€D; for some real J’<A, and satisfying also the weak Neumann 
condition at b*. By the weak Neumann condition at 6* for such a w we mean 
‘or n=2 that uw should possess a continuous extension over Di fe\ |ar| = d}, 
that each first partial of w should be in L, over b<|x|<0-+1, and that 


r—>bt 4 Xj 
|ae|=r lave 


3.9) 0 =lim i w (ae) rail Mx day W040, (x) 

‘or every real-valued function w which is continuous over b6S<|a|<0b+1 and 
Hossesses continuous first partials throughout b<|a|<6-+ 1 which are in L, there. 
[In (3.9), o, denotes »—1-dimensional “‘area’’ measure. For n=1 the weak 
Neumann condition at 0* means simply, with Dj =(— 0, — b) u(b, +0), 
3.10) lim, Wa) = C= Lite dl eg : 

Finally we estimate Nj(A; qo) in (3.7) by the usual separation of variables 
n (3.8) (see [3], Lemma 4.1, plus Addendum). Here for integer »20 we define 


3-11) aly (ie (t) —4(408) Sin 2)? 1 Ae sen — 2) 


over real r=>R, for dimension n= 2, and 


(3.12) Hell eae 
yver such r for dimension »=1. We also define 


(29-+n—2) [y+ n—3 
|» (p=2 


(3.13) LQ =2 At yt, 710) —=1, 
hgh?) = 0 ih 15 LA O\t= 2s 


h,(v) = Mm 2359-6, 
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Define ,N}(A; », g) over real A< 0 and integer y= 0 as the dimension of the closed 
complex linear manifold in L, (b, + 0c) generated by all real-valued w€ L,(b, +) 


which are solutions of 
(3.14) — w"'(r) — pf, (r) w(r) =A’ w(r) 


over (b, +-cc) for some real A’</, the A’ depending on w, and which also have 
a finite limit w’(b*) such that 


(3.45) 0 = w'(b*) — (n — 1) (20) + w(d*). 


With the definitions (3.8)—(3.15) and with g,(x)=—g(|#|), our separation of 
variables statement is that, for }=R,+1 and g as in Condition 2.B, over A<0 


6.10) Ng (Ai ga) =D hn 0) NGA: », 8) 


We indicate as follows the proof of (3.16) cited above (the hypotheses upon 
y= —g are different there than here, but in both cases they imply that the 
real-valued function g is continuous over [b, +o) and in L,(b, +o) and has 
(lim inf {— g(r)})=0, which suffice). For n=2 let {Y,,,}, with 1S m<h, (pv), 

r—-+0o ; S 
denote the L,(S,o,,) orthonormal set of h,,(v) real spherical harmonics in £,, of 
degree v, where S= {9CE,,| || =1} is the £,, unit sphere. Letting J, ,< A, 541<0 
denote the p=1 indexed negative 4’ eigenvalues of (3.14) under (3.15) (each of 
multiplicity 1), and w, ,¢L,(b, +-co) the corresponding normalized real eigen- 
function, we trivially compute that 


(3-47) ty mpl) = [opp (OY, «(lap acy fa ee 


defines u, ,,»CL,(Dj) satisfying (3.8) with 2’=4, , and satisfying the weak Neu- 
mann condition at b* from (3.15). Since the Y,,,, are orthogonal in L,(S,o,) 
for distinct » as well as m, and since (3.14) and (3.15) make the w, , orthonormal 
in L,(b, +c) for each y, the (3.17) functions are orthonormal in L,(D}) for 
m,v, p distinct. Thus the left side of (3.16) is = the right side for n=2. Con- 
versely, letting w be an eigenfunction satisfying (3.8)—(3.10) for 4’, we see that 


(3.18) %,m(t) =r %—D f u(r) Y, ,, (8) do, (®) 


S 


satisfies (3.14) and (3.15) trivially, and thus N=,» and v, ,=cw, , with c= 
(2), m> @y,p) = (4, My mp») Of (3.17). Thus, {Y,, m} being complete orthonormal in 
L,(S,o,), by the Fubini and monotone convergence theorems we have 


»| |u|? du, =f j |u(rD)|? do, (®)} dr 


(3.19) =F 
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where A(A’) = {integer »=0|/, ,=A' for some p} and where in the last line of 
(3.49) p= M with p, for each v€ A(A’), the unique # having 2, ,=2’. Thus (3.19) 
shows that such a w is contained in the closed linear manifold in L, (Dj) generated 
by the functions (3.17), proving equality in (3.16) for n=2. For n=1, (3.16) 
becomes N; (A; qo) =2 Ny (A; 0, g), which is obvious from Dj; = (— 00, — b)u(b, +0) 
with b>0. 

Having established (3.16) for g(a) = —g(|a]|), we use it with (3.7) to obtain 
the boundedness of m(A; Dj, qo) over A<O, which as noted after (3.6) suffices 
to prove our main Theorem 2.C. Here from (2.2) we see in (3.11) that there 
exists some finite integer »,20 in case »=2 such that ,f,(r)<0 over r=R, 
and y>y»,. Thus q(r)=—,/,(r)=0 there shows by application of Lemma 2.E 
to (3.14) and (3.15) with k= (m — 1) (2b) + that ,N; (A; », g) =0 over A<0, b= R +1, 
y>yv,. Hence in case n= 2, (3.16) becomes 


(3.20) Ny (A; Go) = 2» hy, (v) Ne (A; », 8); 


in case n=1, (3.20) with »,=0 is equivalent to (3.16), so (3.20) holds throughout 
A<0 and b= R,+1 with fixed »,=0 for any dimension n. 

Next define ,N,(A; g) exactly like ,Nj(A; y, g) in (3.11) —(3.15), except that 
in (3.14) we substitute g(r) for the ,f,(r) of (3.11) and (3.12). Similarly to (3.1), 
we also define for b= R,+1, ky =(m—1)(20)7 as in (3.15), g=—g or —,f, 


(3.21) Tp (b, 9) = sup int {| (O°) +S (lo (x )|2? + g(x) |w(x)| *) dx} 
with infimums and supremums taken as in (3.1) for D=(b, +e); m,(A; g) is 
the number of integers p21 such that t,(b, —g)<A, and m,(A; ,f,) is the same 
for t,(b, —,/,)SA. Now it is well known ({3], equality (3.27) noted in the 
Addendum last paragraph with here k,=y there) that m,(A;g)=,N,(A;g) and 
My (A; nfy) =rNe (A; », g) over A<O and b=2R,+1. Since (3.11) and (3.12) make 
—_,f,(r)2—g(r) over r=Ry except when »=2 and »=0, clearly m,(A; ,f,)S 
m, (A; g), and thus 


(3222) Ny (A; », 8) SiN, (A; 8) 


over A<0, b= R,+1, integer v=0 except for y=0 when n=2. 


Now 


ris (ms) dssaris*ds=1 “for ree 
shows (2.3)’ implies (2.3), and thus (2.3) always may be assumed over r¢[R,, + 0%) 
for some real R,=>R,+1. Thus g satisfies (2.6) on [R,, +0), and Lemma 2.D 
yields ,N, (A; g) 1 over A<0 and RZ2R, for some real R,=R,. Taking b2=R, 
in (3.22) and (3.20) thus gives our desired Nj (A; 9) SNy<-+ 0 over A<0, pro- 
viding we can also estimate ,N; (A; 0, g) when »=2. 

To handle the exceptional case in which y= 0, n= 2, note that (2.3)’ is assumed 
and of, (r) =2+r-2+-g(r) in (3.11). Hence (2.6)’ is satisfied by 2f), and Lemma 2.D 
yields ,N#(A; 0, g)S<1 over A<O0 and RZR, for some real R;5=Ry+1. Taking 
b=(maxR,, R,), when n=2, thus finishes the proof of Theorem 2.C. 
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4. Discussion of Lemmas 


To prove Lemma 2.E, assume, contrary to the assertion there, the existence 
of a real-valued solution w of (2.8) satisfying (2.9) and 0< af [u(r) }?dr<-+ oo. 


Then (2.8) with g bounded also makes w’’(r)€L,(R, 9); fie it is well known 
([3], equation (3.5) plus references) also that w’(r)cL2(R, co) and that 


(4.1) Fur) ) + [w'(r)]}®} dr = — u(R*) w'(R*) = — k[w(R*))?, 
using (2.9) also. Hence 


0=—kiu( P= J {lw'e (g(r) — A) [w(r)?har 

(4.2) fe = 

= f {a(r) — A} [w(r) Pde S {Lint g(r)] — A} J [w(r) Par > 0, 
R r2R R 


an obvious contradiction. Thus Lemma 2.E follows. 

Turning now to Lemma 2.D, it is well known ([4], we are indebted to Pro- 
fessor HARTMAN for pointing this out to us) that for f as assumed here the con- 
clusion of Lemma 2.D is implied by the equation 


(4.3) — wu" (r) — f(r) u(r) =0 


being non-oscillatory at +-co, by which we mean that there exists some real 
R=R, such that every not identically vanishing real-valued solution wu of (4.3) 
over (R, +00) possesses at most a single zero in (R, +o). But (2.6) is one 
of HILLeE’s sufficient conditions ((5], Theorem 7, Corollary 1 or Theorem 12, case 


p=0) that (4.3) be non-oscillatory at + oo. Likewise, since rfis@ds=}, 
(2.6)’ implies 2 


(4.4) BiG a= arf {> + (sIns)-*}ds 


over r2R,, which is another one ([5], Theorem 12, case P= 1) Obie 
sufficient conditions that (4.3) be non-oscillatory at +-co. Hence Lemma 2.D 
has been proved. 

In closing we remark that other non-oscillation conditions have been given 
for (4.3), but most seem to imply those of HILLE quoted above; naturally each 
such condition leads to a modification of (2.3) or (2.3)’ for which our Theorem 2.C 


follows. If (2.6)’ is dropped and (2.6) is strengthened to fr rf(r)dr<-+os, it is 


fairly easy to give a direct proof (not appealing to non- Eeiaten of (4. 3)) of 
Lemma 2.D; for from the well known asymptotics at oo of solutions of (2.4), 
such Solutions in L,(R, +00) are O(e~®") and satisfy 


co 


(4.5) u(r) eran — f {SRR S— 21) 45) (5) ds, 


from which the violation of (2.5) can be deduced, at least for those HE (0, 2/2) 
given by (3.15). 
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Difference Methods 
for Two-Dimensional Mixed Problems 
for Hyperbolic First Order Systems 


VIDAR THOMEE 


Communicated by LARS HORMANDER 


This paper is concerned with a boundary problem of mixed type fora first- 
order hyperbolic system of partial differential equations in two independent 
variables. 


In two previous papers [10], [11], the author has treated the linear problem 
by L®methods. The existence proof, which is carried out in [//] by use of 
functional-analytic tools, is based on a priori estimates derived in [/0}. However, 
it is not suited for numerical solution of the boundary-value problem. In this 
paper, an alternative existence proof is given by means of an explicit difference 
scheme. This method has the advantages of being applicable for numerical 
purposes and of yielding in a simple way the regularity properties of the solution. 
Moreover, it is suited for the diagonable quasi-linear case. 


The linear problem has been treated previously by other methods by CAmp- 
BELL & ROBINSON [/] and in more than two dimensions by DurF [4]. PEYSER [8] 
has considered the closely connected case of a hyperbolic equation of arbitrary 
order in two dimensions. 

Difference methods have been much used in the theory of partial differential 
equations (see ¢.g. RICHTMYER [9], with references). The success of a difference 
scheme is connected with the possibility of deducing a priori estimates for the 
solutions, or with the “‘stability”’ of the difference scheme. Such estimates are 
also the central tools in this paper. LAx & KELLER [6] have considered an 
alternative difference scheme for a similar boundary problem and proved stability 
under the assumption of the corresponding existence theorem. In a recent paper, 
KreIss [5] has studied the stability problem for a class of equations containing 
the one treated here. 

From the extensive literature on the Cauchy problem for the same system 
the papers by Courant, Isaacson & REEs [2], and by Courant & Lax [3] 
should be mentioned. In [2] the same difference approximations are used as 
here and also are suggested for study of the stability in the case of the mixed 
problem. In [3] a method of existence proof by means of integral equations 
which is used for the Cauchy problem is suggested also for the mixed problem. 
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1. The differential and boundary operators 

Let Y be the set of points {(x, t)|*=0, £20, x+¢/A<1} in Euclidean two- 
dimensional space. Here 4 is a positive constant to be specified later. Let S* 
and S~ be the sets Wn {t=t)} and Wo {x=0}. 

Let (A, B) be the set of mappings from A into B, and let 2(A4, B) and 
€"(A, B) be the subsets of linear and of continuous and q-times continuously 
differentiable mappings. (In case g=0 we usually omit this index.) Let C be 
the set of complex numbers and C, the set of vectors v=(y,..., v,) (in the 
sequal often considered as column vectors) with components in C. Let C,,.,= 
&(C,,, Cy,), .e. the set of complex-valued matrices with , columns and jy rows. 
For M=w or Xz and W=¥, £, S- we put 6h(W)=G1(W, Cm). (We 
often omit M and W and simply 
write @%.) 4 

Assume that 9 is a real dia- 
gonal matrix €@> (VW) and that 
GCC) m(V). The equation 


ou g iu Gu 


(DA) GL = a a 


then defines a linear hyperbolic 
operator LEL(G,(V), Cn (V)). 
We suppose that @ has the form 


D* 0 f? 
Q= ( Fig. 4 


0= D- 


where Y* and Dare diagonal matrices in Gos ymt(V) and Gp-ym-(¥), which 
are positive and negative definite in Y, respectively, and the zeros represent 
blocks with only zeros. Correspondingly we can write u=(u*,u~) where 
ur CEL (V), wu €CEL(V), and partition Y into YEG mr(V) and FEC xm(V). 
With a similar partitioning of Y, (1.1) can be written 


‘ 
Gag oy, 
ot Ox 
ee OU, a Ou - 
eee reaped BY: 


Let 0° and 0 be the restriction operators from Y to /° and -, respectively, 
so that for we Gm (VY) (M=p or fy XM); 


0° u(x) =u(x,0), @ u(t) =4(0,?2). 
On £° we consider the initial-value operator 
Pu=u, LE€2(C(Y), Cn(F)), 
and on S~ the boundary operator 
fu=Bou, CELUCRM(Y), G(F)): 


where & has the partitioning A=(4*, I~) with BS eC as a) and. 1 istthe 
identity matrix of the type m~ xm, With these notations the problem to be 
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considered can be formulated: Given ¥, f°, #. Find uCEm(V) such that 
Gye’. 
(4-2) Piya s®, 
[Uae 
2. The difference operators 


We shall approximate Y by means of a difference operator. We introduce 
for arbitrary 1 the grid 


V ={(Eh, tk)|t=0,... 


0 
Fig. 2. ¢ = points of V; o = points of Vp (n=10; p=2) Fig. 3. ¢ = points of V; o= points of V> (n=10; p=2) 
and its intersections with Y* and SF, 
S'= (Eh, tk) |F=0,...,0— th, OS tS, 
So (OF eB) ir mee lee eee 
We shall need also 
Vi Sith, tk) v=0,..., 081s oA ath 
Vi ={(Eh, tk)|t=0,....n—p; €=0,....n—p — th, 
y= V,o ies = Vea S5; Sp=V,oS*. 
For M=w or yx. and W=V, S° etc. we put @y(W)=%(W, Cw). (We 


often omit the index M.) For ve€@m(V) we now define (for brevity we write 
ve, instead of (Eh, th), etc.) 


—— bi ne 
Agvs = A: Us 7 =h i ie Be vee) 
cs _ f+ 
Ag Ce = A: Urey t 


Here A, and A. EX (€(V) €(V,)) and Az = 
conten ie ue FEX(C(V),@(V,)). The 
definitions of A; and A, for v€G(S°) and @ me ee ce corresponding 


. 
| 
| 
| 
| 
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Let R,, ky,7°, and r~ (pf omitted when =0) be the restriction operators 
mam 1.) -/* and LZ to V,,V,, S°, and S~, respectively, and o® and o~ the 
restriction operators from V to S® and S~ so that for instance if A EG(V), then 
R,AC@(V,) and R, A, ,=A(Eh, th) if (Eh, TR)EV,. 

We shall now define the difference approximations of #* and Y-. We put 
D*—=R,D*, DD=R{D, Gt=R,G*, and G-=R,G- and define for ucG@,,(V) 
f= (urs u-), UCC +(V), u- CE ,-(V)): 


Ltu=A,u*— D* Aput+G*u, 
L-u=A,u- — D- Agu + Gu. 


Here L*EL(G,,(V), Gu (Va) and L-E€L(G,,(V), G,- (V;)). We also introduce 
operators corresponding to f° and f~. We put for u€@,,(V) with B=rB 


(2.1) 


Py =o, LF O=IBG We 


Here €2(G,,(V), G,,(S%)) and -EL(G,(V), Gy-(S>)). 


In addition to the boundary problem (1.2) we shall consider the discrete 
boundary problem: Given F=(F*, F )ES,(V), £°€C, (Sf), £ CCu(S), 
find u€@,(V) such that 

Er Geshe where PreK, F *, 
Ler, where F- = h,F-; 
PU ee Wihere 9/7 LC 


lel if ae WETS) ]t = tee 


(2.2) 


We have the following theorem: 


Theorem 2.1. There is one and only one solution uC@,,(V) of the discrete 
boundary problem (2.2). 


Proof. Theorem 2.1 follows from the following explicit scheme, which is 
easily obtained from the definitions: 


(2.3) Us 9 = 1 Ug, Ea) Mea. Sts 

ogy Mer ADE ates + UP = ADE) the RC ete + AL 

2.4 : d : , 
&€=0,1,....%—T—1; 

0.5) Up ry = (L- + AD; ) U5, — ADE Me_-1,2 — RG Ue, RL Ue 2, 
1. WS Ty 

(2.6) My eg Sl ty — Bh Mo, +41- 


3. Norms 


We define for a vector p=(91,---, P,) €G (W) or a matrix p= (9;;) € Gx u,(W) 
(W=Y, V,, S°, O=O= the point (0, 0), etc.) 


lelw = sup|9;|; lel = sup | ijl - 
1,W 4,J,W 
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Oe 
Let «, and « be non-negative integers, a=(a,,o,), |o|=a,to, and #= 


é i (a) At = Ate Ar. We put 


Ox ot j 
(3.1) ae eee W Bt, Ff? 5 
a 7 a 0 = 
(3.2) l?li,,2= sup |4 P| Wri? W=VeV oe ’ 
a|Sq 


lplog= Sup |@ elo, | Plo, = sup |4* eo: 
 folsa |a| Sa 

In the cases W=£9(S°) and S£-(S-), the sup is taken only over « with a= 0 

and a,=0, respectively. The derivatives in (3.1) are in these cases ordinary 

derivatives. (When p=0 in (3.2), it is usually omitted.) 

We shall state some simply proved properties of these norms. 

Lemma 3.1. For each q, there is a constant C,, independent of n, p we of 
a,€€,(W), 1=1, 2; W=V,, Ve, Sp, Sp 
such that 


——F 


| a as|W,¢ Ss C,|4|w,¢ |aelw,¢- 


2 


/ 


SSN ORES Lemma 3.2. If a€¢"(W), W=7, 
ENENSS V,S, and SF, respectively, then 
SENSES 
[Sie SSS |R, a lv, sa la\ra 
SIS 7 | 
PESTON OSSIS DRESS TESS [Rs aly-, a= |4|x93 
IS ES SEIS 0 Z 
ESSE SSeS) Se ot hss «= Ary 

Fig. 4 


7p @|s- @S la|¢-¢: 
For the functions A€@y(V,) we define an extension operator EE L(y (V,), 


€u(V)) in the following way. We first extend (if necessary) A to the whole 
of V by putting 


A, ,=Ao np for t= Pa oy es Sis iS eaceg hls == bie 


Then we triangulate V by means of the 3(m-+1) straight lines x/h=j; t/k=]; 
x[h-+-t/kR=7;7=0, ..., n. In each of the triangles the value of EA is determined 
by the values in the three edges so that EA(Eh, tk)=A,, and so that EA is 
linear in the triangle. Similarly, we define the extension operator E in the case 
A-€€y(V;). In this case we extend A~ to the whole of V by putting 


Ag, =A, p12, for t=0,....m0—-p—1, §=n—p—t41,....n—7, 


hp = AG opts for UND. teas We Ci 0 ore om ates 


and then proceed as before. In the same way we define extension operators 
E°ER (Gy (Sp), Gu (F)) and E-EL (€m (Sp), Gm (S-)) by linear interpolation such 
that E°a°(£h)=a? and E-a-(rk)=az where the definitions of az and az are 
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extended through 


a, =A»; tT=n—p+l,...,N. 


Some of the important properties of the extension operators are described 
in the following lemmas. 


Lemma 3.3. If A€Gu(V,), A-E@u (Vp), @€ 6 (S2), and a-~CGy (Sz), then 
JEAly=|Aly,; [EA |y=|4y; |E°@|po=|@|s0; | E-a|y-=|a"|s:. 
Lemma 3.4. If a€ @,(W), then for n> 00, we have in $y(W), W=VV, SF, S, 


EAR, a> 0a; © la) =9; 
EA* Roa > a; lal <4; 
E° A* Pa > a; itl = Ge p= 0; 
E-A*ra—> a; lala Wag=s0: 


(Convergence in @y(W) means uniform convergence in W.) 

Lemma 3.5. If a€@j(V), then |ERa—aly<4(h+hk)? |alys. 

Lemma 3.6. If a; ,,C€@(V,), 1=1, 2, n=1, 2,... and Ea; ,>a;in @(V) when 
n—> oo, then E(a, ay ,) >a, 4, in 6 (V) when n->oo. The corresponding statements 
hold for V,;, S$, and S;. 

The following two lemmas will be used in the existence proofs. 


Lemma 3.7. I} for every n, u,€Em(V) and |u,|y4<C, where C is a constant, 
then a sequence {n,\{° can be found such that Eu,, converges to a function u in 
€Cu(V) when 7-0. 

Lemma 3.8. If for every n, u,€@m(V) and |uny|y,qi1SC, where g>0 and C 
is a constant, then a sequence {n,}1, can be found such that Eu,, converges to a 
function u in 6y,(V) and E A*u,, converges uniformly to o*u for |x| <q when] >. 

Lemmas 3.7 and 3.8 are consequences of ARZELA’s theorem (cf. PETROwsKI [7], 
p.2784f.). 

4. Estimates 
In Sections 4—6 we shall make the following assumption about the constant A: 


0<AS|D|/. 


The main result in this section is the following theorem. 


Theorem 4.1. There is a constant C, independent of n and u€G,,(V), such 


that 
|uly SC{|Lt uly, + [Lo uly + [2a] se + [bu] sz}. 


For the proof we need Lemma 4.2, which is proved by aid of Lemma 4.1. 
Lemma 4.1. Assume that a,=0, T=0, ..., n, and 


G44 4a,+271C,(@,+ kK), +=0,1,....%—1. 
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Then for 6>0 and tSOn, 
aata,t (rj 


Proof. We have 
Gray t KS (t+ OQ70 G+ 5) se" (a+), 
and hence for t<0n, 
g, + Kae! (ag pK). 


Lemma 4.2. Assume that at>0, a;=0, a,=max(aj, az), T=0,...,n and 
that 
at.,sattntc,(a4,+4), ang UAE Ft yee 


1 max{a; +n71C,(a,+ K); C3(at.,+K)}, = 0,1, 2 ake 


at 
1 


(4.1) 


Then there exists a constant C,, independent of n, such that 
(4.2) a, = C,(4) + K), T= 0400: 


Proof. Without restricting generality we may suppose that 


(4.3) C, 2 C,/Cs. 
We put 
b6 =a, 
by max tD. 502), TC Ome ier 
(4.4) oe : 
pe Oe nt Ce), SSK bier ek, 


br4. = max{b; +220, (b,+.K), C(t +K)}, 1=0,1,...,.0—1. 


It follows at once that af<bF, a,Sb,. Since b)=ap, it follows that (4.2) is 
proved by the following inequality: 


bee C (BAL Ry) aera 


Suppose that for some t in (4.4), b;=C,(b}+K). Then we also have 
b741=C3 (6341+ 4), for according to (4.3) we have for this t 


6, +n2+*C,(b,+ K) = CA + n(C,/C,) (0, + K) + kK} 
= C,{b? +n1C, (b,+ K) | i C3 (bt41+ K). 
It follows that there is a t) such that 
b41= 6, ++4C,(0,+ K), To ty: 
6741 = C3(07,1+ 4), TZ=%. 


The case where b, ++ C;(b; + K) for r=0,..., also fits in here if we agree to 
put t)= 7 in this case. For t< 1, we thus have 


bo Sb, 4 UCH(b OR) 


According to Lemma 4.1, since ty <n, we have 


(4.5) b, Seb 4K) Se (0, HW CoG i vee ey 
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We obtain 
(4.6) bt 41 = df, +04, (0, +.K) Cy (b,, + K). 


For t >t) we have 
Ory Sb; +0 7C, (6; + 0; + K) SOF +04 (1+ C,) (6 + 4), 


and according to Lemma 4.1 
brs Cobia + A); 

that is, because of (4.6), 

OFS ONO IS) | E Tee 
Since 

b =C, (0 + kK), > Tes 
it follows that 

STO (0 aap x0 a oe 


and hence according to (4.5) 
b,<C(b+ 4K), PSE 
which completes the proof of Lemma 4.2. 
Proof of Theorem 4.1. If for w€@,(V) we put 
Ge =|e' lox, a= |e |sz, @,=|elsc, 


K=|Ltuly,+|Louly- + |b uls-, 


then from (2.3), (2.4), and (2.5), with appropriate C,,C,,C, we obtain the 
inequalities (4.1), and thus the theorem follows from Lemma 4.2. 


Theorem 4.1 implies the following theorem for the differential operator. 
Theorem 4.2. There is a constant C, independent of uCG)(V), such that 
(4.7) |ulyS CLL uly + |lou| go + |u| y-}. 
Proof. We have according to Theorem 4.1 
|Ru|yp SC{|L*Ruly, + |L-Rulp- + | Rew go + | Ru|s-}. 
Therefore because of Lemma 3.3 
JERulySC{EL*Ruly+|EL Ruly+|E°PRu| go +|ET Ruly-}, 


and hence (4.7) follows from Lemma 3.4 and Lemma 3.6. 


For later use we state the following two lemmas, which are obtained by 
examining the proofs of Lemma 4.2 and Theorem 4.2. 


Lemma 4.3. Assume that at =>0, a7 =0, a,=max (az, a7), T=0, ..., n, and that 


Gp a, nC (ao), T= 0; 4 tS 


(iy max {a; +n1C,(a,+ 5); C3(at41+ K,)}, == 0) 41,5..,0 45 
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Then there exists a constant C,, independent of n, such that for OS 01 
at <C,(a+6K,+ 6K), 


Lemma 4.4. Let o<6<1 and V8=VO(0StS 64}, ( =Sf-AV*. Then 
there is a constant C, independent of uCG,(V), such that 


| e*]y3 SC{|Ou| go + d|Lulyst olf u|.y-y}- 
|u| S Cf [Lm] po + 8] L ul y9 + [ful gr}. 
Theorem 4.3. Assume that G,G,BEE!. Then there is a constant C, inde- 
pendent of n and of uEG,,(V), such that 
(4.8) [uly SC{|L uly g +L 4 |y-g + [Pes ¢ +12 u|s-¢ +|#lo,ct- 


Proof. We shall prove the theorem by induction over g. For g=0, the 
statement is that of Theorem 4.1. One easily proves the following representation 
of the successive differences with respect to T: 


At uf 1 =A De 74; Ab vig, + -ADE r+j) 4) ud +k Pee oO; 
(4.9) AL We ool (LD, om) A Ug yoAD, ip Ale at a i kiy a E 
A! 7 Uo, eae (1 uy. +1 — Bi41 40, Brey 


Here Pos ,=L* ue ,—GE Ue ,, and pane 


Amey re Abed ap aa Dew AL (DE = (4 UE A GE ete — Us, ake 
We have 
(4.10) | Bis,2| SC {Ata «| + | eva + [Le aly, aie 


Now assume that the theorem is proved for g=7;—1. (4.9) and (4.10) then give 
Ae | set S [AL a] sp + C (AE wl g¢ + [Lt ely, ; + IL-l yz,i), 

| Ay | spn S max {Ab a] gx + Cy (| At | se + |L* uly, j 

Cy (|A% u*| sea + [o> a 5-, + |\L*u 


mi)? 


ae |L-u 


Bele 

and thus according to Lemma 4.2 

(4.11) |At uly, C{\Atu set |Lt uly i+ [L-ulp- 5+ |o-4|s-,)}- 

Now using (2.1), we easily obtain 

(4.12) [Ap m|sp SS C{[P a] 0,¢ + [L aly, ja (Lo ml ye ga + [4 o,j + | ly, ;-a}- 


Using successively the formulas (2.1) and, in the neighborhood of the line 
—+t=n, the easily proved formula 


A gue, = (I~ +A Deva > 1) Ag uz 1tAd,ug , — 
**: AGE 41, t—1 Met 7-1 oh AL” Ue ya EE Os, 


(4.13) 


emneaiedl 


————— 
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we obtain for 0</<7 
(4.14) |APAP ul, SC{|Ai aly + [Lt ely, pa + [Lo | ye pa + | ly ja}. 


The induction assumption (4.11), (4.12), and (4.14) now prove (4.8). 
Theorem 4.3 implies the following theorem for the differential operators. 


Theorem 4.4. Assume that 9,G, BCE. Then there is a constant C, inde- 
pendent of uC Git1(W), such that 


(4.15) [Hyg SC{| Luly, + [ful yo t+ |f-uly-,}- 


Theorem 4.4 is proved in the same way as Theorem 4.2. The term |u| , 
which is obtained from the term |u|) , in (4.8) can be estimated by means of 
the terms on the right-hand side of (4.15). 


5. Compatibility conditions 
We say that FE@ (VY), PEGS), £-CE(S-) satisfy compatibility 
conditions of order q at the origin (for the boundary problem (1.2)) if there 
is a uCGi(V) such that 


(5.1) a (Lu) (0, 0) = & F(0, 0), lal Sq—1, 
(5.2) # (0x) (0) =# £0), BS, (4=Z), 
(5.3) #C-#) (0)=# 4-0), BS4q, (4=—,). 


If for instance F, #°, and # vanish in the neighborhood of the origin, then they 
satisfy compatibility conditions of order g for every g at the origin. We have 
only to put #=0 in ¥.. In order to define the compatibility conditions it would 
only be necessary to require regularity of F, #°, and # at the origin. However, 
this fact will not be emphasized since the global regularity will be needed in 
the existence proofs. 


We obtain at once 


Theorem 5.1. A necessary condition for the solvability of the boundary problem 
(1.2) in Gi(V) ts that F, £°, and ~~ satisfy compatibility conditions of order q 
at the origin. 

In the same way as in the proof of the Cauchy-Kowalewsky theorem we 
realize that for arbitrary FY and /° we can find a polynomial # such that (5.1) 
and (5.2) are satisfied. The equations (5.1) and (5.2), however, determine 
6*a(0, 0) for |x| <q, and thus d°f-(0), Bq, in (5.3) cannot be chosen arbitrarily. 
The following theorem is easily proved. 


Theorem 5.2. There are matrices M,, and N,,, depending only on D, Y, and 
B, such that if F, £°, and ~~ satisfy compatibility conditions of order gq at the origin, 
then for ySq 


(5.4) d’ £-(0 =2Mp d® £°(0)+ +2 ~oF (0,0), 


78 VIDAR THOMEE: 


and such that if on the other hand (5.4) is satisfied for F, £°, and ~~, then we can 
find a polynomial %=1 (x, t) of order q such that (5.1), (5.2), and (5.3) are satisfied. 


In particular, for g=1 we have the following theorem. 


Theorem 5.3. Assume that BCE him-(S). Necessary conditions for the 
solvability in G1(¥) of the boundary problem (1.2) are the following relations : 


(5.5) (0) = A(0) £9(0), 
(5.6) &f-(0) = B(0)[A(0, 0) £9 (0) — (0, 0) £9 (0) + F(0, 0)] + 4B(0) £°(0). 
Proof. We have according to (1.2) 
~ (0) = B(0) u(0, 0) = BO) £°(0), 
df-(0) =4 (Blt) 40, )) 0 = B(0) $* (0, 0) + €B(0) «(0, 0) 
— B(0) P00, 0) £% (0,0) — ¥(0, 0) #(0, 0) + Luo, 0)| + dB(0) u(0, 0) 
= B(0) [B(0, 0) d£° (0) — ¥(0, 0) £9 (0) + F(0, 0)] + dB(0) £9(0). 


It turns out that (5.5) and (5.6) are also sufficient conditions for the solvability 
of the boundary problem (1.2). This will be proved later by means of Lemma 5.1. 
In the existence proof we consider the discrete boundary problem (2.2). 


Lemma 5.1. Assume that 9, G9, FCO; B, $°, £ €C* and that the conditions 
(5.5) and (5.6) are satisfied. Then there 1s a constant C such that if u,, is the solution 
of the discrete boundary problem (2.2), then for all n 


[tploa=C. 
Proof. We have according to (2.1) 
(5.7) A,ug,= Di, Ai ui,— Gi Ue + Liu, ., (6,7) = (6,0),.141,0) Men ie 
In particular, we obtain 
(5.8) Avis = Dia kG eer ee 


and hence for = 0 owing to the regularity assumptions we obtain the boundedness 
of A,u¢,o. Because of (2.6) for r=1, 2 we have 


According to (5.5), (5.9) is valid also for t=O. Therefore we get 

(5.10) A, Uo,0 = A, fo — A, (Btu) 9. 

Because of the regularity assumptions, Lemma 3.2 and the already proved 
boundedness of A, 1), the boundedness of A, Uo, follows from (5.10). From 
(5.8) we get , 

(5.11) A; A, 46,0 = A; (D* A; f9*)4 — Ag (G* f)\y + A; Foro 
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and hence the boundedness of A;4,u¢,.9. From (5.7) we obtain 
A? U6,0 = A,(D* A; U*)o 9 aaa A,(B* u)o 4 ae A, Fo'o, 


and hence the boundedness of A? 19 9 follows from (5.11), the boundedness of 
A,uo,o, and the regularity assumptions. From (5.9) we get 


At Uo,0 = Az fo — Az (Bt u*)o 9, 
and hence the boundedness of A? Uo 0 follows as above. 
It remains to prove the boundedness of A;A,u9,9. We have 
A, 43,0 = Dio 4 Uy 0 — Gio, 9 +L my 9 =DioAsfo —Gioht+ Fo. 
Together with (5.8) and (5.10), this gives after a simple computation 
A; A, Uo, 9 = n{B(0) [B(0, 0) d# (0) — F(0, 0) £°(0) + F(0, 0)] + 
+ dB(0) £°(0) —df-(0)$+O(1) when noo. 


According to (5.6) this implies the boundedness of A;4,u9 9, and the proof of 
Lemma 5.1 is completed. 


Lemma 5.2. If FEGI(V), PEG (S), EECA (S-) and if 


OF (0, 0)= 0, Nees es 
(5.12) “zee (0) eee =O B99 
a’ £-(0) =0, BSq, 


then there 1s a constant C such that af u,, is the solution of the discrete boundary 
problem (2.2), then for all n 


(5.13) logit cs C- 
Proof. According to (2.3), (2.4), (2.5), and (2.6) we can write 
la| |o| 
a Croat 19s bausig aRiG + nlel- ‘ > R, _ Ff, Et? 


E4+7S|o|—1 
where the coefficients P:, Q,, and R,, are bounded when noe. According to 
(5.12) we therefore have 
AU =O (niat i Ol) ewheny 7-00 for |x| = 94-1, 
which is (5.13). 
Remark. It might be expected that in general if FE@{(V), LC Gf" * (F), 
£-€i*'(S-) satisfy compatibility conditions of order g, then (5.13) would be 


valid. For g=1, this was the case according to Lemma 5.1. For g>1 this is 
not always the case, as is proved by the following example. 


Consider the boundary problem (m=1) 
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The compatibility conditions of order 2 are here 


g(0) =/(0), 
S(O — 0), 
£5 (0) = f(s 


A simple computation for the solution w,, of the corresponding discrete boundary 
problem gives 


A? A, %, 9 = 2" (f'(0) + g’(0)) +n(Ag'(0) — f’’(0)) + O(1) when ”—>0©o, 
A, A? up, 9 = 0? (f'(0) + g’(0)) + 2 ((A — 2) 80) +f/(0)) +O(1) when n—>oo. 
It follows that if A<1, then |#,|) 3=O(1) when m—co only if f’’(0)=g"(0) =0. 


6. Existence and uniqueness theorems 
We first state the uniqueness theorem. 
Theorem 6.1. If F, £°, £ €@, then there 1s at most one solution uc G2} (Vv) 
of the boundary problem (1.2). 
Proof. This follows at once from Theorem 4.2. 
We now prove the existence theorem in @!(¥) and first consider the special 


case q=1. 


Theorem 6.2. Assume that 9,9, 2B, F, f°, £- CC? and that the conditions (5.5) 
and (5.6) are satisfied. Then there is a solution uC €,(V) of the boundary problem 
(e2). 

Proof. Consider the discrete boundary problem (2.2). Owing to Lemma 3.2 


we have 
Pela 


Vy,2 We Fig ste t fesse C 


Lemma 5.1 gives for the solution wu, of (2.2), |wnlo,2<<C. Therefore according 
to Theorem 4.3 we have |w,|y,.<C. Using Lemma 3.8, we can find a sequence 
{m,jr such that Eu,, converges uniformly to w€G,(Y) when k->co and such 
that EAgu,,, EAgu,,, and EA,u,, converge uniformly to du*t/éx, du-/@x, and 
du/ot, respectively. We assert that w is a solution of (1.2). Lemma 3.4 and 
Lemma 3.6 give 


wu = lim EL 4, =m ER, Ft aF*, 
R—>0o k—>0o 
£-u=lm EL-u, =lim ER} F- =F, 
k—oo k—> oo 
onset Oy eet a 
£94 = lim E Pu, = lim E° 7° Psa £0, 
tu = hm Eo lu, = lim. Borg =f. 
which proves Theorem 6.2. 


Theorem 6.3. Assume that q=1 and QF, G, B, F, f°, ¢£- €6""" and that F, 
f°, f satisfy compatibility conditions of order q at the origin. Then there is a 
solution uC; (>) of the boundary problem (1.2) 
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Proof. We first assume that the conditions (5.12) are satisfied. Consider the 
discrete boundary problem (2.2). The assumptions and Lemma 3.2 give 


ete Fr lyejqgta + | /°|s°, 41 a [fl s- gta =. 
| Lemma 5.2 gives for the solution u, Of (2.2) 

| %n lo, ¢+1 SC. 

} According to Theorem 4.3 we therefore have 

ltnly gi S C. 


Using Lemma 3.8, we can find a sequence {n,\{° such that Ew,, converges uni- 
formly to uc@i(Y) when k->oo and such that EA%w,, converges uniformly 
to o*u when k->oo for |a|<g. The function w is as in Theorem 6.2 a solution 
of (1.2). 

We now assume that the conditions (5.12) are not satisfied. Since F, £°, f- 
fulfill compatibility conditions of order g at the origin, owing to Theorem 5.2 
we can find 7€ G2"? (VY) such that (5.4), (5.2), and (5.3) are satisfied. We con- 
sider the boundary problem 


Lipa — Soy vines, 
Eo Jo bos = ODE S.0, 
(EF A on TS 


Here the conditions (5.12) are fulfilled, and thus there exists a solution ue G/(V ). 
The function w+ u€G?(V) is then a solution of the original boundary problem 
(1.2). This proves Theorem 6.3. 


Theorem 6.4. Assume that D9, G9, B, F, f°, £- ©", and that the conditions (5.5) 
and (5.6) are satisfied. Then we have for the solution u, of the discrete boundary 
problem (2.2) 

E A“u, > &u, lic diy 


where u is the uniquely determined solution of the boundary problem (1.2). 


Proof. We prove the assertion for |«|=0. For |«|=1 it can be proved 
analogously. It follows as in the proof of Theorem 6.2 that every sequence {n,}%° 
contains a subsequence {n,}{° such that Eu,, converges to a solution up 0 of 
the boundary problem (1.2) when k->oo. According to the uniqueness theorem 
(Theorem 6.1) we have Upmpjo=u, where u is the only solution of the boundary 
problem. Hence we have Eu, —>u when noo. 

It is also possible to obtain results concerning the rate of convergence, for 
instance the following theorem. 

Theorem 6.5. Assume that, G, B, F, £°, $- €E*, and that compatibility rela- 
tions of order 2 are satisfied at the origin. If u and u, are the solutions of the 
problems (1.2) and (2.2), then 


zt -1 
|Eu,—uly~=O(n4) when noo. 
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Proof. According to Theorem 6.3, ucG2(v). We can write 


|Eu, —UlyS|E(t,— Rwy +|ERu— uly =|, — Ruly + |ERu— uly. 


According to Lemma 3.5 we have |ERu —u|y=O(n*) when n— oo. To estimate 
|x,, —Ru|y we use Theorem 4.1 and obtain 


|u,—Rulps Cla, ae |L-u, —L~Ru|y-+ 
+ |l°u, —1°Ru|so+|l-u, —I- Ru|s-} 
= C{|Ltu,—L*Ruly, +|L-4,— L-Rul|y-}, 


since u, =lRu=/? and I-u,=l-Ru=f-. We have (0S0;=1, + =1, 2) 


\L+u,¢7—L*Ru, | = |F (EA, th) —L* Ry, ,| 


out out + + 4 Avy, + iG* 
| Gg + Gru —A,ut++QDt Atut—Y aaharn 6 


aut out ee} 
= | ey dl S Agape & 


SE (Eh, (¢-+0) 8) — DEA, th) SE (E48) h, Th) 
SF (R+|D yh) |w'|y2=0(n) when n—>oo, 


and similarly for the minus case. This proves the theorem. 


7. The quasi-linear case 
In this section we consider the quasi-linear boundary problem 


Ou Ou “ oh ' 
a g~. =F, Q=Y(x,t,u), F=F(x,t,u) in V, 
(7.1) Ua LO. Oils eee 


Usa P"\ OMe VOM, hn a Pea). 


Vv, S*,and S-, with A to be specified shortly, have the same meaning as before. 
About 9, F, #, - we make the following assumptions: 


1. £°CG,(S). There are closed sets (u=(uy, ..., tm), U= (ty, ---, Um+)) 
2={e1wl0S¢et Oa f— a. |u — £9(x)| < et}, 
Q={Ew)|0Stsd,, |u*— gor0)| <e} 
such that FC Cy (Q), DEC hx m(Q), £- €G2(Q-). In Q, D can be written 
Salles 


where D*>0 and D-<0. 
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2. The following compatibility conditions are satisfied at the origin: 
a) £-(0, £°*(0)) =£°-(0), 
B) “Zp (0,42*(0)) +-Fe= (0,210) |2°(0,0, £2(0)) “4 (0) +-F*(0, 0, 24(0) 
= D(0, 0, £9(0)) of (0) + F-(0, 0, £°(0)), 


| where Of |du* is the matrix with columns Of /0%;,, KAA, ..., mM. 
Under these assumptions we shall give a local existence proof for the boundary 


problem (7.1). Let A be a positive number <|QJ|j' and let h, k, V, etc. have the 
| same meaning as before. We introduce 


YEE NOSt= Oh, Vi ae 
We consider the discrete boundary problem 


Aru — Di Arun, =Le,;, on Vz 

Alu, — De, Aru. = Te on Vi 
Wage |e sO S° 
Ue he) se Ola So. 


| Here Dry QE, th, u(Eh, tk)), etc. The following result may be obtained as 
in the proof of Theorem 2.1. 


Theorem 7.1. There is one and only one solution u€€,,(V) of the discrete 
boundary problem (7.2). 
The solution can be defined by induction over t by means of the following 
scheme: 
Weis fe, E=0 
Ue ADE Mey et ADE) ui, RE, . €=0,..,n—t—1, 
Ue roy = (I + ADS) Ug — ADE Mg_-1,2 +h Fe, cea 


een 


(7.3) 


yreeyN—Tt—1, 
Uo, 741 = peat 
Theorem 7.2. Assume that 9, F, ~°, £~ satisfy 1. and 2. Then tf u ts the 
solution of the discrete boundary problem (7.2), there is a d>0 and a constant C 
such that independently of n, |u|yo 9<C. 
For the proof we need Lemma 7.2, which is proved by means of Lemma 7.1. 


Lemma 7.1. Assume that a,=0, t=0,1,...,”, and 
G44 4,+n"f(a,), B= 0}; coop = WW, 


where f(x) is non-decreasing, continuous, and >0 for x20. Then there exista d>0 
and a continuous, non-decreasing function w(x) defined for x= 0 such that a, Sy (a) 
for rn 1S. 

6* 
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Proof. We put b,=a)+211f(a). We have b)=4 and 
(7.4) bevy — by — mf (b,) = 7 [2 (dg) — f (do + 2f (ao) 7 2)]. 


Consider the continuous function x (2)=2f (a) —f(a+2f (49) z). We have 
1(0)=f (a) >0, and therefore there exists a 6>0 such that y(z)=0 for all z 
with 0<z<6. With this 6 we have according to (7.4) 


(7.5) bx Zd, +mtf(h), tH S6. 


Since f(x) is non-decreasing, (7.5) implies a,<6,, and therefore with p(x)= 


x+26 f(x), 


A, Sa) +206 f (a) = y (4%), TUS 0, 
Lemma 7.2. Assume that at =0, a7=0, 4,=max (aj, az), T=0,..., 0, and that 
Api G, 1 0 *7(G,), T=0,1,...,.%—1, 


az 41 max{a; + nf (a,), g(at,1)}; tO tc2yG ee 


where f(x) 1s non-decreasing, continuous and >0 for x= 0 and g(x) 1s non-decreasing 
and continuously differentiable with g'(x)=1 for x20. Then there exist a d>0 
and a continuous, non-decreasing function p(x) defined for x=0 such that a,<y (ap) 
[00.7 WSO. 


Proof. We put 


b= aj 
b= max (O25 ba), Cie (eee ees 
(7.6) 
Oia 0, see 7(0.),, 8 Oils Sesto == ihe 


by 41 = max {b; + nf (b,), g( tia}, t=0,4) [jae 


Since f and g are non-decreasing it follows easily that a*<60#, a:b. Tie 
therefore sufficient to prove that there exist a 6>0 and a continuous, non- 
decreasing function p(x) such that b,<y(by) for tn4<6. Suppose that for 
some t in (7.6), b;=g(bt). Then we also have b7,,;=g(b%,1) for (0<8S1) 


8 (r+1) = 8 (br + m7 (by) = (bt) +n F(b,) g' (ds + On f(b) = be + nf). 
It follows that there is a ty such that 


bia = 0; cata | (0) aren. 
(7.7) Or 41 = 8 (br 41), TZ. 


The case where by + g(b7) for r=0,..., also fits in here if we agree to put 
T) =n in this case. For t< ty we thus have 


bi Sb, + nf (b,). 
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According to Lemma 7.1 there exist a 6,>0 and a continuous, non-decreasing 
function y,(x) such that 


(7.8) b, Sy, (b), TNO, tao. 
We have 
(7.9) Be 41 = 3, +07 F(b,,) Sd, + fra) S Yr Oo) +f (v1 (b0)) = Ye (60) - 
For t >t) we have 
bry = 0; +017 (b,) SbF + n+ f (max (8}, g(63))). 


The function f(max(x, g(x))) is non-decreasing, and therefore according to 
Lemma 7.1 and (7.9) there exist a 6, > 0 and continuous, non-decreasing functions 
ws (x) and y,4(x) such that 


(7.10) Op Ws (07, 4.1) = Ws (Wo (bo) = 4 (00); TUESSO55 T= Tor 


Finally, (7.7), (7.8), and (7.10) imply with 6=min (6, 6,) and y(x«) = max (y,(x), 


pa(*), & (pa(x))) that 
b, Sy (%), Tite 20, 


which completes the proof of Lemma 7.2. 


Proof of Theorem 7.2. From the difference scheme (7.3) one easily deduces 
the following representation of the successive differences with respect to T: 


At ue 41 = ADE 45 Ab Mba + (It — ADE 4) Aue +R Ofer, 
f— Ole ee, 4 F 7 — 1S 
(7.11) At Ug eet = (ae ADs, 1+) Ae Us — ADs 245 4 Mga, +h Qj, &,2> 
fF=1,...,.n—t—J—1; 
Ae Uo, 141 = AD ee 
Here 05. ,—F,, and 
Ofis,62 =A, Ofte + 4, Di 14; Ai {(DE 2) 1 (4,4. — FE}. 


In the scheme (7.3) we consider only &, t such that (£h, tk, ug ,) EQ that is, 
—,r such that O0<rk<6,; OSEhS1—Th; |u,—fe|Se. It is then easy to 
prove that 


(7.12) |Q¢2-| GC, (|4,u|& +1), 
} | Ose] SC (| a|s7 +1), j>4, 


where C; is bounded with [At «| 5, 1=0,...,7—1. (In fact we shall not use 
these estimates for 7 >2.) Therefore from (7.11) and (7.12) we find that 


|A, u*| sr S|A, u*|57 + RC, ([4, u[$ +1), 
|A, #-|gr1S max {|A, u-|o¢ + RC, (|4, u|$e +1), Ky (lA, u*|s¢+1 + 1)}, 
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where the assumption that the constant Ky satisfies K,=1 is no restriction. 
We shall see now that |4,u|s., 1e., A,“éo is bounded. For the + case and 
£=0,1,... and for the — case and =1, 2,... this follows at once from (7.2) 
and assumption 7. Further we have in view of the compatibility condition 2a 


A, Uo, 0 = A, fo; 
which is bounded with A,1,9. Lemma 7.2 now shows that there exist a constant 
C and a 6,> 0, 6,=06,, such that 
|A,v|or SC, TES 0s. 


The assumption |v; , —/?| Se for r<t%, where t9k< do, thus implies |A,u|szSC 
for t<ty). It is therefore an easy induction to prove that 


lu e—A[SC-h<e, tk d,=min(Z, d). 


From (7.2) and (4.13) we can estimate 4.u,,, and we thus find that there 
exist a 6 >0andaC such that |u | yois C. It remains to estimate the second-order 
differences. We consider A? ue and first observe that Aug 6 is bounded. For 
A2uzy,E=0,1,...and A?ug 9, &=2,3,... this follows from (7.2). Also A?up > =A? fo 
is bounded with Aius,o, 4=0,1,2. Finally using (7.2), we find that A? Ux 9 is 
bounded with 4;4,u9 9. We have easily from the definitions 


A, A, %,0 = 1 9-(0, 0, £°(0, 0)) Liane (0) + F-(0, 0, £°(0)) — es (0, £°*(0)) — 


oF” (0, £2*(0)) [a(0, 0, £2 (0)) 24 (0) + F*(0,0, £° (0))|}4+-0 (4) when n—>0o , 


out ae 


and thus according to the compatibility condition 26, A:4,u99=O(1) when 
n—>co, We have with a constant K, 


|A? u*|sr1 S |A? w*| 2 + kC,(|A? «| 52 oh. 1) 
[Az w|sz1 S max {|APu-|gr + & Cy (|Az ule +1), Ky (|? u*|sz1 + 1)}, 


and therefore, since |A?u|,. is bounded according to Lemma 4.2, |A?u|y,8 is 
bounded. In the same way as in Section 4 we establish the boundedness of 
|Azu|y,3 and |A,A,u|,,s, and thus Theorem 7.2 is proved. 

As a consequence of Theorem 7.2 we have the following theorem, which is 
proved in the same way as Theorem 6.2. 


Theorem 7.3. Assume that D, F, ¢°,é satisfy 1. and 2. Then there exists a 
6>0 such that the boundary problem (7.1) has a solution in V°. 


We supplement the existence theorem by means of the following uniqueness 
theorem. 


Theorem 7.4. Assume that £°€60(S°), that 


Q = {(x,t,u)|OS*<1, 0O<tS4, | uw — £9(x)| <8}, 
Q={(u')|0StS 4, [ut — £*(0)| <e}, 
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and that F €G;,(Q), BEE. (Q), LC SL(Q-), and 


les re 0 

Oma 
in 2, where Dt>0O and D-<0. Assume further that u and v are solutions in 
€n(V°*) of the boundary problem (7.1). Then there exists a 6>0 such that u=v 


in V°. 
Proof. Let 6,>0 be so small that (x,t, w(x, t))€Q and (x, #, v(x, t))€Q 
for all (x, t)CV. Let Lw= 2? — Qu) ce and f-w=w-(w=(w*, w-)). Then 
w=u —v satisfies the linear boundary problem (¥ (u) =F (x, t, u(x, t)), £-(ut)= 
f(t, u*(0, 4))) 
Lw=Fu) — Fr) + (Bw) — Bo) 2, 
w=, 
fw =u) — £0). 

According to Lemma 4.4 we have for 6< 63 


jo|y0s C{| F(x) — Fo) + [Bu — BON] IF 0) — Fw). 


and thus with 
0D 
Ou 


Ov 
Ox 


ag- 


out 


| 
I 
Q 


i + 
ys | 


? 


Q 


OF 
=m 
O | Ou OF 


where the definitions of the matrices 04%/du and 0Y/du are obvious, we get 
(7.13) |w*|yo< CC, d|wlys, 6 < 63. 
Further, in the same way for 6< 6, Lemma 4.4 yields 
|w|ys SCC, dl wlyo+CC,|w*|y, 
and thus using (7.13), we have 
| w-|y-e< (CC, + (CC)?) dl wlys, O=aOn 


end with C,=—CC, -+ (CG)? 
|wlyos C,6|w|ys. 


Now if 6<min(6,, Cz"), we obtain |w|y-e=0, which proves the theorem. 
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